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ON CLASSICAL SOLUTIONS FOR VISCOUS POLYTROPIC FLUIDS 
WITH DEGENERATE VISCOSITIES AND VACUUM 

YACHUN LI, RONGHUA PAN, AND SHENGGUO ZHU 


Abstract. In this paper, we consider the three-dimensional isentropic Navier-Stokes 
equations for compressible fluids allowing initial vacuum when viscosities depend on 
density in a superlinear power law. We introduce the notion of regular solutions and prove 
the local-in-time well-posedness of solutions with arbitrarily large initial data and vacuum 
in this class, which is a long-standing open problem due to the very high degeneracy 
caused by vacuum. Moreover, for certain classes of initial data with local vacuum, we 
show that the regular solution that we obtained will break down in finite time, no matter 
how small and smooth the initial data are. 


1. Introduction 


In this paper, we investigate the local-in-time well-posedness and formation of singu¬ 
larities of classical solutions to the compressible Navier-Stokes equations for isentropic 
flows when viscosity coefficients, shear and bulk, are both degenerate and the initial data 
are arbitrarily large with possible vacuum states. The system of compressible isentropic 
Navier-Stokes equations in R 3 reads as 

Pt + di v(pu) = 0, 

( pu)t + di v(pu ®tt) + VP = divT. 

We look for local classical solution with initial data 

(p,u)\ t=0 = (p 0 (x),uo(x)), iGI 3 , 

and far field behavior 

(p, u) —>• (p, 0) as |a;| —>- +oo, t > 0, 

where p > 0 is a constant. 

In system (11.111 . x = (xi,X 2 ,xs) € R 3 , t > 0 are space and time variables, respectively. 

p is the density, and u = (t/ 1 ', vP\ € R 3 is the velocity of the fluid. We assume 

that the pressure P satisfies 

P = Ap\ 7 > 1, (1.4) 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 
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where A is a positive constant, 7 is the adiabatic exponent. T denotes the viscosity stress 
tensor with the form 

T = p(p)(Vu + (Vu) T ) + A (p) divii I 3 , (1.5) 

where I 3 is the 3x3 identity matrix. We assume in this paper that 

p{p) = otp 5 , A (p) = (3p 5 , (1.6) 

p{p) is the shear viscosity coefficient, A (p) + | p(p) is the bulk viscosity coefficient, a and 
/3 are both constants satisfying the following physical constraints 

a > 0, 2a + 3/3 > 0. (1.7) 

Furthermore, in this paper, we require that the constant 5 satisfies 

1 < 5 < min { ^ ^ , 3 j. (1.8) 

In the theory of gas dynamics, when we derive the compressible Navier-Stokes equations 
from the Boltzmann equations through the Chapman-Enskog expansion to the second 
order, cf. Chapman-Cowling |7j and Li-Qin [23], we find that the viscosity coefficients are 
not constants but functions of absolute temperature in a power law. For isentropic flow, 
this dependence on the absolute temperature, from the laws of Boyle and Gay-Lussac, can 
be reduced to the dependence on the density, that is to say, the viscosity coefficients are 
proportional to powers of density, see [23] , 

When inf-rpo^) > 0; it is well-known that the local existence of classical solutions 
for (ll.ll) - (jl.3H can be obtained by a standard Banach fixed point argument due to the 
contraction property of the solution operators of the linearized problem, c.f. Nash [32] . 
However, this approach does not work when vnl x p${x) = 0, which occurs when some 
physical requirements are imposed, such as finite total initial mass, finite total initial 
energy, or vacuum appearing locally in some open sets. 

When viscosity coefficients p and A are both constants, for the existence of solutions 
of the three-dimensional isentropic flows with generic data and vacuum, the main break¬ 
through is due to Lions [23], where he established the global existence of weak solutions 
provided that 7 > |. This result is improved to the case 7 > | by Feireisl-Novotny- 
Petzeltova[13j. and even to non-isentropic flows in [13, [15]. However, the regularities of 
these weak solutions are fairly low, and the uniqueness problem is widely open. On the 
other hand, the local well-posedness problem in higher regularity class with possible vac¬ 
uum initial data requires comprehensive treatment on the degeneracy in time evolution 
in momentum equations (I1.1IB . Since the leading coefficient of ut is p which vanishes at 
vacuum, there are infinitely many ways to define velocity (if it exists) when vacuum ap¬ 
pears. Mathematically, this degeneracy brings forth an essential difficulty in determining 
the velocity when vacuum occurs, since it is difficult to find a reasonable way to extend 
the definition of velocity into vacuum region. Physically, it is not clear how to define the 
fluid velocity when there is no fluid at vacuum. A reasonable framework was proposed by 
Cho-Choe-Kim [ 8 j ID, Ilf through a compatibility condition 

- divT 0 + VP(po) = \/P09 (1-9) 

for some g £ L 2 . In turn, a local theory of strong or classical solutions with initial vacuum 
was established successfully for three-dimensional case; see also Duan-Zheng-Luo [12] and 



NAVIER-STOKES EQUATIONS 


3 


Luo |29j for two-dimensional case. More recently, Huang-Li-Xin }20j extended this local 
classical solution an]) to a global one with small energy and vacuum for isentropic flow 
in M 3 . In these results, the uniform ellipticity of viscosity plays a key role in the a prior 
estimates of higher order terms of velocity u. 

When viscosity coefficients p and A are both density-dependent, system m has re¬ 
ceived extensive attentions in recent years. Instead of the uniform elliptic structure, the 
viscosity degenerates at vacuum, which raises the difficulty of the problem to another 
level. A remarkable discovery of a new mathematical entropy function was made by 
Bresch-Desjardins [2J for A (p) and p(p) satisfying the relation 

A(p) = 2(n'(p)p- n(p)). (1.10) 

For example, in our problem, such an entropy structure (11.1011 exists when /3 = 2a{5— 1) > 
0. This new entropy offers a nice estimate 

p'(p)YleL°°([0,T]-,L 2 (R d )) 

VP 

provided that //(p 0 )^= £ L 2 (R d ) for any d > 1. This observation plays an important 
role in the development of the global existence of weak solutions with vacuum for system 
m and related models, see Bresh-Dejardins El H, Mellet-Vassuer ED, and Vassuer- 
Yu [ID, and some other interesting results, c.f. [5, 0, 221 E3 [28, 39j SOj. However, we 
remark that, in spite of these significant achievements mentioned above, the local well- 
posedness of classical solutions in multi-dimensions with vacuum is still open due to the 
high degeneracy of this system near vacuum region. Our existence result in this paper is 
a solid step toward this direction. 

From (11.4|1 - (I1.5I) . we notice that 

divT = —p s Lu + V// • §(u), 

where the so-called Lame operator L and the operator § are given by 

Lu = —aAu — (a + /3)Vdi vu, §(u) = a(Vu + (Vu) T ) + /3 divul 3 . (1.11) 

When constructing the classical solutions with arbitrarily large data and vacuum, in ad¬ 
dition to the issue shown above for the constant viscosity case on the degeneracy in time 
evolution in momentum equations, there are still two new significant difficulties due to the 
appearance of vacuum. The first one lies in the strong degeneracy of the elliptic operator 
divT caused by vacuum. We note that under assumptions (II.61) - (1 1 .8f) . viscosity coefficients 
vanish as density function connects to vacuum continuously in some open sets. This de¬ 
generacy gives rise to some difficulties in our analysis because of the less regularizing effect 
of the viscosity on the solutions, and is one of the major obstacles preventing us from uti¬ 
lizing a similar remedy proposed by Cho et. al. mmm for the case of constant viscosity 
coefficients. The second one lies in the extra nonlinearity for the variable coefficients in 
divT due to (11.61) . We emphasize here that, unlike the constant viscosity case, the elliptic 
part divT is not always a good term in the regularity analysis for the higher order terms 
of the velocity. For example, if we want to get the estimate on || V 3 k,||x / °°([o ,t] ; l 2 (k 3 )) inde¬ 
pendent of the lower bound of the initial density, we need to deal with an extra nonlinear 
term 


div(9^/9 5 (a(Vu + (Vu) T ) + /^divula)) 
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where ( = (Ci, ( 2 , C3) is a multi-index with |£| = C1+C2+C3 = 3- Therefore, much attention 
will be paid to control this strong nonlinearity. 

In order to overcome these difficulties, some new ideas have to be introduced. In [27], we 
obtained the existence of the unique local classical solutions for system with vacuum 
far field under the assumption 


<5 = 1, ck > 0, a + ft > 0 


( 1 . 12 ) 


for two-dimensional space, aiming at the application to shallow water models. In [27] we 
observed that, assuming p > 0 , the momentum equations can be rewritten as 

u t + u-Vu + +Lu= (¥£) -S(«). (1.13) 

7 — 1 V p J 

Then, (11.131) implies that if we can control the special source term ( 77 ) • § (u) when 

vacuum appears, the velocity u of the fluid can be governed by a strong parabolic system. 
However, this result only allows vacuum at the far field. The corresponding problem with 
vacuum appearing in some open sets, or even at a single point is still unsolved. 

Motivated by this observation in when p > 0 , for our case under the assumption 
m, instead of (11.131) . (ll.ljl o can be rewritten as 


ut + u ■ Vu + p s 1 Lu 


2A~/ 

5-1 


c p 


27 —6 — 1 s-l 

g ~ 1 V p 2 + 


25 6-1 6-1 

j—jp 2 Vp 2 


S(u). 


(1.14) 


Then if we pass to the limit as p —>• 0 on both sides of (11.141) , we formally have 


ut + u ■ Vu = 0 when p = 0. (1.15) 

So (jl . 14D - (j 1.15 f) imply that actually the velocity u can be governed by a nonlinear degen¬ 
erate parabolic system when vacuum appears in some open sets or at the far field, which 
is essentially different from the parabolic system in (11.131) in the sense of mathematical 
structure. Based on this observation, we introduce a proper class of solutions and prove 
the local-in-time well-posedness of solutions with arbitrarily large data and vacuum in 
this class for system (11.11) using a new approach which bridges the parabolic system (11.141) 
when p > 0, and the hyperbolic system (11.151) when p = 0. 

In order to present our results clearly, we first introduce the following definition of 
regular solutions to (ll.ll) - (ll.3D . 


Definition 1.1 (Regular solutions). Let T > 0 be a finite constant. ( p,u)(x,t ) is 
called a regular solution in M 3 x [0, T] to the Cauchy problem 11.111 - 11.3\) if 

(T) (p, u) satisfies the Cauchy problem 11.1\) - 11.3\) in the sense of distribution; 

(B) p> 0, p^-^eC([0,T];ff 3 ), (pX) t eC([0,r];H 2 ); 

(C) u E C([0, T];H S ') D L°°([0, T]; H 3 ), p^r V 4 u € L 2 ([0, T]; L 2 ), 
UiGC([0,T];R 1 )nL 2 ([0,T];D 2 ); 

(D ) ut + u ■ Vu = 0 holds when p(t, x) = 0, 

for any constant s' € [2,3). 
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Here and throughout this paper, we adopt the following simplified notations, most of 
them are for the standard homogeneous and inhomogeneous Sobolev spaces: 

I/Ip = II/IIlp(r 3), ll/IU = ll/IU s (R3), l/b = ll/llo = II/IIl 2 (r3), 

D k ' r = {/ € L} oc (R 3 ) : |V fe /|r < + 00 }, D k = D k ’ 2 , \f\ Dk , r = || f\\ D k,r^ 3) (k > 2), 

D 1 = {/ £ L e (R 3 ) : |V/| 2 < 00}, | f\ D i = ||/|| d i (R 3 ). 

A detailed study of homogeneous Sobolev spaces can be found in m- 

Remark 1.1. This notion of regular solution for compressible Navier-Stokes equations 
was first introduced by Yang-Zhu m for one-dimensional case. We defined the regular 
solutions in m for the case 6 = 1. Compared with m, a significant difference is that the 
admissible initial data in this paper are much broader. Actually in m, in order to make 

sure that the source term ■ §(w) appearing in the equations hl.lfft) is well defined in 

H 2 space, we need 

K 2 ), ^e^nh'nfl^R 2 ), 

p 

which means that the vacuum must and only appear in the far field. In this paper, we only 
need 



So the vacuum can appear in any open set or in the far field. 

Now we give the main existence result of this paper. 

Theorem 1.1. If the initial data (pq,uq) satisfy the following regularity conditions: 

Po > 0, ^p 0 2 - p L ^,u 0 '\ € H 3 (R 3 ), (1.16) 

then there exists a positive time T* and a unique regular solution ( p,u)(x,t) in R 3 x [0, T*] 
to Cauchy problem (Cup-(Cnp. 

Remark 1.2. The weak smoothing effect of the velocity u in positive time t € [r, Tf, 
Vt e (0, T*), tells us that the regidar solution obtained in Theorem \l.l\ is indeed a classical 
one in (0, T*] x R 3 . The details can be found in the Appendix. 

Remark 1.3. We can also consider the flow with external force in the momentum equa¬ 
tions ( 11 . 1 b ) such as 

(pu)t + div(pu <g) u) + VP = divT + ph. (1.17) 

Assume that 

h € C([0,T];P 1 (R 3 ))nL 2 ([0,T];P 3 (R 3 )) (1.18) 

for some T > 0, then the same conclusion as in Theorem, 11.11 still holds if we replace the 
condition ( C) in Definition 11.11 with 

ut + u- Vu = h when p = 0. 

To achieve this result, one only needs to make minor modifications on our proof for The¬ 
orem 1.1. We will not go into details on this matter. 
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As a direct consequence of Theorem 11.11 and by the standard theory of quasi-linear hy¬ 
perbolic equations, we have the following additional regularities for some specially choices 
of power 6. 

Corollary 1.1. Let 1 < 5 < min j|, j, or l 5 = 2 (7 > 3). If the initial data (po,uo) 
satisfy II. 161) . then the regular solution ( p , u ) obtained in Theorem. 11.11 also satisfies 

p-peC([0,T*);tf 3 ), pt £ C([0, T*); H 2 ). (1.19) 

Compared with previous results for 6 = o!m and 6 = imi our analysis in the proof of 
Theorem 1 1.1 1 is quite different. From the observations shown in (I1.14I) - (I1.15I) . the behavior 
of the velocity u near vacuum under the assumption (| 1 , 8 [> is more hyperbolic than the 
cases in mi and m ■ This phenomenon leads us to developing a different approach in 
the analysis on the regularity of u. As a matter of fact, in B uni HD E 3 , the uniform 
ellipticity of the Lame operator L defined in (11.111) plays an essential role. Therefore, one 
can use standard elliptic theory to estimate \u\ D k+2 (0 < k < 2 ) by the 19'-norm of all 
other terms in momentum equations (I1.1IL for the case 5 = 0 in [S HfliCD], or in (11.131) for 
the case 5 = 1 in [271 [42j. However, for our case, it is clear from (11.141) that the standard 
parabolic theory only offers the following weighted estimate: 

^^6L 2 ([0,T];L 2 ), for |C| = 1,2,3,4. 

This estimate alone is not enough for the regularity analysis on u and p. In order to obtain 
the desired estimate 

u € L°°([0, T], H 3 ) n C([0, T]; H s ) for any s' E [2,3) 
shown in Definition D we turn to the help of symmetric hyperbolic structure 

ut + u- \7u. 

These estimates are accomplished along with a vanishing viscosity limit argument. 

Another challenging question in the theory of fluid dynamics is the problem of global 
regularity. This problem is extremely difficult for three-dimensional compressible Navier- 
Stokes equations, especially when the initial density contains vacuum. When viscosity 
coefficients are constants, global regularity was achieved when initial energy is small, see 
m for the case away from vacuum, and m for the case allowing possible initial vacuum 
for the solutions in the class of [ 8 ]. In particular, the results of |20| showed that no 
singularity will form in finite time for smooth solutions in the class of P if the initial 
energy is small. On the other hand, there are a lot of finite time blowup results obtained, 
say 0 29,133], [37, 38] SQ], for various solution classes or conditions. However, it is not 
clear yet if such solutions exist locally in time. For our case, when the viscosity coefficients 
satisfy the conditions (ll. 6 D - (ll. 8 l) . we will show that the solution we obtained will break 
down in finite time for certain classes of initial data with local vacuum, no matter how 
small the initial energy is. This is in sharp contrast to the case of constant viscosity as 
shown in [sum m- This is because our problem behaves more closely to compressible 
Euler equations in vacuum region, due to strong degeneracy of the viscosity. 

In fact, we will present two scenarios for singularity development in finite time from 
local vacuum initial data. The first one is motivated by Xin-Yan [38], called the initial 
data with isolated mass group. 
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Definition 1.2 (Isolated mass group). (po(x),uo(x)) is said to have an isolated mass 
group ( Ao,Bq), if there are two smooth, bounded and connected open sets Aq C R 3 and 
Bq C R' 3 satisfying 

' a 0 CB 0 C B Ro C R 3 ; 


* p 0 (x ) = 0, VieB 0 \4 



po(x)dx > 0; 


( 1 . 20 ) 


for some positive constant Rq and constant vector uq S R 3 , where B Rq is the ball centered 
at the origin with radius Rq. 


Our first blowup result shows that an isolated mass group in initial data guarantees the 
finite time singularity formation of regular solutions. 

Theorem 1.2 (Blow-up by isolated mass group). If the initial data (po,uo)(x) 
have an isolated mass group (Aq,Bq), then the regular solution ( p,u)(x,t) in R 3 x [0, T m \ 
obtained in Theorem 11.11 with maximal existence time T m blows up in finite time, i.e., 
T m < Too. 


For the second scenario, we explore the hyperbolic structure for the system in vacuum 
region. For this purpose, we introduce the following concept. 


Definition 1.3 (Hyperbolic singularity set). The non-empty open set V C R 3 is called 
a hyperbolic singularity set of (po, uq)(x), ifV satisfies 


j p 0 (x) = 0, VieV; 

(sp(Vuo)nR-4 0, V®€ V, 

where Sp(\7uo(x)) denotes the spectrum of the matrix Vuq(x). 


( 1 . 21 ) 


This definition is inspired by the global existence theory of classical solutions to the 
compressible Euler equations in [T 8 j and [Mj, for initial data satisfying the following 
conditions 

• po is small and compactly supported, 

• Sp(Vuo) n R _ =0, V a; € R 3 . 

Our second blowup result confirms that hyperbolic singularity set does generate singu¬ 
larity from local regular solutions in finite time. 


Theorem 1.3 (Blow-up by hyperbolic singularity set). If the initial data (po,uo)(x) 
have a non-empty hyperbolic singularity set V, then the regular solution ( p,u)(x,t ) on 
R 3 x [0, T m ] obtained in Theorem 11.11 with maximal existence time T m blows up in finite 
time, i.e., T m < Too. 


We now outline the organization of the rest of this paper. In Section 2, we list some 
important lemmas that will be used frequently in our proof. In Section 3, we first refor¬ 
mulate our problem into a simpler form in terms of some new variables, and then we give 
the proof of the local existence of strong solutions to this reformulated problem, which 
is achieved in four steps: ( 1 ) we construct global approximate solutions for a specially 
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designed linearized problem with an artificial viscosity r/ 2 Lu in momentum equations; (2) 
we establish the a priori estimates independent of artificial viscosity coefficient ry; (3) we 
then pass to the limit ry —>• 0 to recover the solution of this linearized problem allowing 
degeneracy in the elliptic operator appearing in momentum equations; (4) we prove the 
unique solvability of the reformulated problem through a standard iteration process. Sec¬ 
tion 4 is devoted to proving the finite time blowup shown in Theorems ll.2til.3l Section 5 
is the appendix where we prove the fact that the regular solution we obtained in Theorem 
II . II is indeed a classical one in (0, T*]. 

Finally, we remark that our framework in this paper is applicable to other physical 
dimensions, say 1 and 2, with some minor modifications. This is clear from the analysis 
carried out in the following sections. 

2. Preliminaries 

In this section, we list some important lemmas that will be frequently used in our proofs. 
The first one is the well-known Gagliardo-Nirenberg inequality. 

Lemma 2.1. [21 j For p £ [2,6], q € (l,+oo), and r € (3,+oo), there exists some generic 
constant C > 0 that may depend on q and r such that for 

f£H l ( R 3 ) and g £ L q (R 3 ) n D 1 ’ r (M 3 ), 

we have 

\f\ p P <c\f\t p)/2 \^f\?~ 6)/2 , 

( 2 . 1 ) 

Iffloo < C , |5|5 (r - 3)/(3r+9(r - 3)) |V 5 |f /(3? ' + ' ?(? '- 3)) . 

Some commonly used versions of this inequality are 

\u\e<C\u\ D l, Moo < C\u\l |Vtt|| , \u\oo < C\\u\\ w l,r. ( 2 - 2 ) 

The second lemma is some compactness results obtained via the Aubin-Lions Lemma. 

Lemma 2.2. [35] Let Xq C X C X\ be three Banach spaces. Suppose that Xq is compactly 
embedded in X and that X is continuously embedded in X\. Then the following statements 
hold: 

(1) IfG is bounded in L p ( 0, T ; Xq ) for 1 < p < +oo, and is bounded in L 1 ( 0, T\X\), 
then G is relatively compact in L p (0,T] X). 

(2) IfG is bounded in L°°(0,T;Xo) and ^ is bounded in L p (0,T; X\) for p > 1, then 
G is relatively compact in (7(0, T; X). 

The following four lemmas contain some Sobolev inequalities on the product estimates, 
the interpolation estimates, the composite function estimates, etc., which can be found in 
many literatures, say Majda m- We omit the proofs of them. 

Lemma 2.3. [30] Let constants r, a and b satisfy the relation 

- = —f - , for 1 < a, b, r < +oo. 
r a b 

Vs > 1, if /, g € W s ’ a n VF s,b (M 3 ), then 

IV s (/ ff ) - fV s g\ r < c s (|v/| a |v s -M| b + |V7Ws|a), 


(2.3) 
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|V s (/ 5 ) - fX s g\ r < C a (|V/| a |V s -^| 6 + \X s f\ a \g\ b ), (2.4) 

where C s > 0 is a constant depending only on s, and X s f (s > 1) stands for the set of all 
partial derivatives dif of order |£| = s. 

Lemma 2.4. [3DJ If functions /, g G H s and s > then fg G H s , and there exists a 
constant C s depending only on s such that 

\\fg\\s<C s \\f\\ s \\g\\s. 


Lemma 2.5. [30] If u G H s , then for any s G [0, s], there exists a constant C s depending 
only on s such that 

i_y l_ 

Hulls' < C'sII^IIq s IMI/■ 


Lemma 2.6. [30] 

(1) If f, g G H s n L°° and |£| < s, then there exists a constant C s depending only on 
s such that 


\\d C x (fg)\\ s < C s (\fUV s g\ 2 + b|oo|V s /| 2 ). (2.5) 


(2) Let u(x) be a continuous function taking its values in some open set G such that 
u G H s (~l L°°, and g(u) be a smooth vector-valued function on G. Then for s > 1, 
there exists a constant C s depending only on s such that 


V s g(u) 1 2 < 


dg_ 

du 


\u 

s — 1 


s —1 
oo 


V s n| 2 . 


( 2 . 6 ) 


The following lemma is a useful tool to improve weak convergence to strong convergence. 


Lemma 2.7. [30] If the function sequence {rc n }^ =1 converges weakly to w in a Hilbert 
space X, then it converges strongly to w in X if and only if 

IHIx > lim SUPn^cc\\w n || X- 


3. Existence of regular solutions 


In this section, we aim at proving Theorem 11.11 and Corollarv ll.il To this end, we first 
reformulated our main problem ( 1 1 . 1 1 ) - (1 1.3 1) into a more convenient form. Since we will 
repeat the integration over M 3 for many times, we will adopt to the following simplified 
notation 

[ f= [ fdx, 

J J R3 

throughout this paper without further specification. All other integrals will be specified 
when they appear. 

3-1 _ 

3.1. Reformulation. Setting <j> = p 2 , system (11.11) can be rewritten as 


ft + u ■ VQ + - - 1 0 di vu = 0 , 

2 Ay 2r — 5—l 9 

ut + u ■ Xu + A 1 


(3.1) 


6-1 


+ f A Lu = Xf 2 • Q(u), 


where 


Q(u) 


<5 — 1 


( a(Xu + (Vu) T ) + /Jdivnla) = - ^ S(m). 


(3.2) 
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The initial data are given by 

{(j),u)\t=o = (0o,uo)(x) = (^p 0 2 (x),u 0 (x)j , x E R 3 , 
with the far field behavior 

(0,u) -» (4>,0) = (p“,0) , as \x\ -A- Too, t > 0. 


(3.3) 


(3 - 4) 

To prove Theorem ll.il we first establish the following existence result for the reformu¬ 
lated problem (13.1 1) - (13.41) . and then in subsection 3.6 we will show that this result indeed 
implies Theorem ll.il 


Theorem 3.1. If the initial data (0o,uo)(x) satisfy 

00 > 0, (0o ~ 0, u 0 ) € H 3 , (3.5) 

then there exists a positive time T* and a unique strong solution (0,u) on R 3 x [0, T*] to 
Cauchy problem h3.1\) - [3.4\) , that is, (0,u) is a solution of the Cauchy problem h3.1\) - [3.4\) 
in the sense of distribution and satisfies 

0-0£ C([0,T*]; LT 3 ), 0 t E C([0, T*]; H 2 ), 

u E (7([0, T*]; H s ) PI L°°([0, T]; H 3 ), <)V 4 «ei 2 ([0,T,];l 2 ), (3.6) 

u t E C'([0,T*];ll 1 )nL 2 ([0,T*];T> 2 ), 
for any constant s' E [2,3). 

We will prove this theorem in the subsequent four subsections 3.2-3.5 as explained in 
the introduction. 

3.2. Linearization with an artificial viscosity. In order to construct the local strong 
solutions for the nonlinear problem (|3.1I) - (I3.4I) . we first consider the following linearized 
problem: 

X ^ _ 2 

0t + V • \70 H- 0&YVV = 0, 


A.'W 2-y—2 o o o 

u t + v • Vu H-V</> 5 ~ 1 + ((fr + rf)Lu = V(fr • Q(v), 

7-1 

(0,u)\t=o = (0o(x),uo(x)), X E R 3 , 

k (0, u) —>• (0, 0), as |x| —^ Too, t > 0, 

where E (0,1) is a constant and 

6 S 

Q( v ) —-(a(Vr + (Vv) T ) + /3divul 3 ) = -—-§(u). 


(3.7) 


(3.8) 


Here 0 and v = (t/ 1 ), v^ 2 \ t/ 3 ))" 1- E R 3 are given functions satisfying the initial assumption 
(0,v)(t = 0, x) = (<^> 0 ,^ 0 ) and the following regularities: 

0-0eC([O,T];H 3 ), 0 t eC([O,T];H 2 ), 

v E C([0,T]; H s> ) n L°°([0, T]-,H 3 ), 0V 4 v E L 2 ([0, T]; L 2 ), 

v t € C([0, T]; H 1 ) n T 2 ([0, T]; D 2 ), 


(3.9) 
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for any constant s' € [2, 3) and fixed time T > 0. 

Such a linearization is specially designed to serve our purpose. Unlike E3, here we 
have to keep the positive and symmetric hyperbolic operator in the momentum equations. 
More precisely, we linearize 


ut+u- Vu -A ut + v ■ Vu, 

instead of ut + u ■ S/u -A ut + v ■ Vu as in m, which plays a key role in the analysis on 
the regularities of velocity u. Generally speaking, in the following proof of the existence, 
we will make a full use of the transport property of the hyperbolic part ut + v ■ Vu. 
On the other hand, the strong diffusion term ( q i 2 + q 2 )Lu helps to ensure the global 
existence of the linearized problem. Moreover, we remark that there exists a carefully 
chosen compatibility between the last term divu on the left hand side of (I3.7D i and 
the viscosity term (</> 2 + q 2 )Lu in (13.71) 9. Because of the lack of positive lower bound for 
<f> o, it is not visible to show that there exists some constant C independent of q such that 


u llz, 2 ([0,T];U 4 ) < C. 


Therefore, in order to obtain a finite upper bound of |^>|_c >3 independent of q from the 
continuity equation (13.7|) i. we need a good estimate like (fiV 4 v € L 2 ([0, T]; L 2 ). For our 
solution ( (j),u ), ||</>V 4 u|| L 2 ([ 0 iT ]; L 2 ) < C can be obtained by the smoothing property of the 
artificial viscosity term (0 2 + r] 2 )Lu. It is not clear if other linearization approaches will 
achieve the same goal, this one (13.71) succeeded. The details can be found in Lemmas 

E2K1 

As mentioned above, a rather standard method gives the following global existence of 
a strong solution (</>, u) to (13.71) for each fixed q > 0. 

Lemma 3.1. For each fixed rj > 0 and any T > 0, assume that the initial data (4>o,uo) 
satisfy \3. 51) . Then there exists a unique strong solution (4>,u)(x,t) in M 3 x [0,T] to \3. 7| 1 ? 
that is, a solution satisfying the Cauchy problem B> with the regularities 

fi-fi€C([0,T};H 5 ), (j) t £ C([0, T]; H 2 ), 

u G C([0, T]; H s ') D L°°([0, T]; H 3 ), cjN\ E L 2 ([0, T]; L 2 ), (3.10) 

u t GG([0,T];i7 1 )nL 2 ([0,T];T» 2 ). 

Proof. The first equation of (13.71) is a linear transport equation, the existence and regularity 
of a unique solution <f> in M 3 x [0, T] to the first equation of (13.71) can be obtained by the 
standard theory of transport equation. 

After knowing <f> from (|3.7|) i . for q > 0, (13.71) 9 becomes a linear parabolic system of 
u. Then it is not difficult to solve u by the standard parabolic theory. Here we omit the 
details. □ 

In the following two subsections, we first establish the uniform estimates independent 
of q, then we pass to the limit q -A 0. 
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3.3. Uniform a priori estimates independent of i]. The main task of this subsection 
is to establish some local (in time) a priori estimates for the solution (</>, u ) to (|3.7D obtained 
in Lemma [3.II independent of the artificial viscosity coefficient r/. For this purpose, we fix 
a T > 0 and choose a positive constant Co large enough such that 

2 + + |<^o|oo + ||0o — 0||3 + ll^olb < Co- (3-11) 


We now assume that there exist some time T* E (0, T) and constants q (i = 1,2,3,4) 
such that 

1 < c 0 < Cl < c 2 < c 3 < c 4 , 


and 


sup (IV'lL + H(t) - <j)\\l + ||v(t)||i) <cf, 

0 <t<T* 

rT* 

sup + \v{t )\ 2 D 2 + \v t (t)\l) + / (\ipV 3 v\ 2 2 + \v t \ 2 D i)dt <cl, 

0 <t<T* JO ' 7 

ess sup (|V^(*)lx>i + K*)li)3 + |ut(i)li)i) + [ (|'0V 4 u|| + |u t |^ 2 )d; t <cj, 

0 <t<T* Jo V 7 


sup <cl 

0 <t<T* 


(3.12) 


We shall determine T* and c* (* = 1,2,3,4) later, see (13.641) below, so that they depend 
only on Co and the fixed constants p, a , /3, 7 , A, 6 and T. 

Let (c f>,u)(x,t ) be the unique strong solution to (13.71) in R 3 x [0,T], We start from the 
estimates for 0. Hereinafter, we use C > 1 to denote a generic positive constant depending 
only on fixed constants p, a, /3, 7 , A, 6 and T, which may vary each time when it appears. 


Lemma 3.2. Let (< f>,u)(x,t ) be the unique strong solution to fl 3. 7| ) on R 3 x [0,T]. Then 

1 + (j) + |0(^)|^o + II0(^) — 0ll§ <C* c §, 

< Cc t \Mt)\ 2 D i < Cci \Mt)\ 2 D 2 <Cc% 
for 0 < t < T\ = min(T*, (1 + cz)^ 2 )- 
Proof. Let f> = 4> — <j). Then from (|3.7I) 1 . we have 

4>t + v • V(ft + ^ ^ 0di vv = 0. (3.13) 

Applying the operator di (0 < |C| < 3) to ()3.131) . we obtain 

{d^(j))t + v ■ = — (d£(v ■ Vcjb) — v ■ Vd|</>) — ^ ^ d^^divv). (3.14) 

Multiplying both sides of (|3.14p by dif>, and integrating over R 3 , we get 

~l^ll < c (|divu|oo|d^|| + A 1 101 2 + A||^0| 2 ) , (3.15) 

where 

A i = I diiy ■ V0) - v ■ VdH | 2 , A 2 = |^(0divv)| 2 . 
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From Lemma 12.31 letting r = b = 2, a = +oo, f = v, g = V0 in (12.311 . we obtain 

Ai < C c (|Vw|oo|Vl c l0| 2 + \VtUdivU) . 

Now we consider the term A^. When |£| < 2, it is easy to show that 

A^<C c (h/>|oo + ||VV>|| 2 )IM| 3 ; 

when |£| = 3, we have 

As <C (|'0V 4 t’| 2 + |V0 ■ V 3 u | 2 + |V 2 0 • V 2 u| 2 + |V 3 0 • Vv| 2 ) 

(|- 0 V 4 ^| 2 + liv^lbll^lls) ■ 

Combining (13.15D - (13. 18D . we arrive at 

^II0(*)I| 3 <c (jM| 3 || 0 || 3 + (Moo + IIV' 0 || 2 )IMI| 3 + IV’V 4 ^^) 

<c(c 3 || 0||3 + c 2 + |0V 4 0 2 ) . 

From Gronwall’s inequality, we have 

||0(0|| 3 <(|| 0 o - 0|| 3 + c|t + J |0V 4 u| 2 ds^ exp(Cc 3 t). 
Therefore, observing that 

J | 0 V 4 u(s)| 2 ds < | 0 V 4 u(s)||ds^ <Ccztz, 

we get 

||0(f) - 0 || 3 < Cco for 0 < f < Ti = min(T*, (1 + c 3 ) -2 ). 

The estimate for 0j follows from the relation: 

4> t = —v • V0 ———0divu. 

For 0 < f < Ti, the following estimates hold 

|0t(f )| 2 <C'(|u(f)| 6 |V0(f )| 3 + 10(f)loo|divv(f)| 2 ) 

<Cc 2 , 

|0t(f)lz?i <C(|u(f)| oo |V 2 0(f )| 2 + |Vu(f)| 6 |V0(f )| 3 + |0(f)|ooV 2 u(f)1 2 

+ |Vu(f)| 6 |V0(f)| 3 ) 

<Cc 2 , 

|0t(f)M <C'(0(f)| oo |V 3 0(f )| 2 + |Vu(f)| oo |V 2 0(f )| 2 + |V 2 u(f)| 6 |V0(f )| 3 

+ |0(0|oo|V 3 u(f)| 2 + |V0(t)|oo|V 2 u(f)| 2 + |V 2 0(f)| 2 |Vu(f)| oo ) 

<Cc 2 3 . 

These complete the proof. 


(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


(3.21) 


□ 
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Now we turn to the estimates for the velocity u in the following two lemmas. In view 
of (11.81) . we define 

4y — 4 

V = <5-1 “ 8 ' 

The following lemma gives some lower order estimates for the velocity u. 

Lemma 3.3. Let (cj),u)(x,t) be the unique strong solution to 3. 7\) on R 3 x [0, T\. Then 


\u(t)\\l<Ccl, \u(t )\ 2 D 2 + \u t (t)\l + J (j</>V 3 u(s)|| + \u t (s)\ 2 Dl ^j ds < Cc 


- - cl K 


for 0 < t < T 2 = min(T*, (1 + C3) 2K ). 

Proof. The proof will be carried out in two steps. 

Step 1. We estimate \diu \2 when |£| < 2. Taking di on equation (13.7f) o . we have 


(diu) t + v • V(<9 iu) + (c 1 i 2 + rf)Ldiu 

= - + v </> 2 • d^Qiy) - (d^fv -Vu)-v V(d<u)) (3.22) 

- (dx((</> 2 + V 2 )Lu ) - {cj? + ri 2 )Ld$.uj + (^(V^ 2 • Q(v)) - V</> 2 • d$.Q(v)j . 


Multiplying (13.221) by d^u on both sides and integrating over R 3 by parts, we have 


1 d 


n~n\ d i u \2 + a\^ (f 2 + r] 2 V (diu)\l + (a + /3)\ y/cf 2 + r/ 2 div^u|: 


2 dt 


= ~ J(v -V(d£u)) ■ d£u - j (y{(j) 2 + rj 2 ) • Q(d£ix)) • d$u 
A \ J -diu + J Vcj) 2 • diQ{v) • cf x u 


7 - 


+ 

7 


J • Vu) - v • V(^u)) • d Q x u 

/ ■ AM)-v^-*«.))•*. 




(3.23) 


i= 1 


Now we estimate the last part of (13.231) term by term. Using Holder’s inequality, Lemma 
IQ and Young’s inequality, we have 
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h = ~ J (vV^u))-^ 

<C'|V?;| 00 |9£u|! < Cc 3 \d$.u\l, 

h = - J (v(</> 2 + T] 2 ) • Q(a£u)) • 

< C 7 | 0 V (^«)| 2 | V 0 | oo |^«| 2 
<C| A /^+^V(^)| 2 |V0| oo |^| 2 
<ylW 2 + ?? 2 V(5 |u)|| + Ccgl^wll. 

For the term 1%, it is obvious that 

/ 3 =- J Vc>|0 ■s- 1 • d^u = J d^cj) s- 1 divc^tt. 

Since 1 < <5 < yy-, > 4. Then for |£| = 0, one has 

2/4'V ( 2-y-«-l 

4=- —Y j 

<mj^ iv0| 2 |u| 2 

47 - 2 , 5-2 

<c , |0| oo 5 - 1 IV^II + CH 2 
<Cc^ +C\u\l 

When 1 < |£| < 2, we have 

h =^—y J di<j)~s=^diwd^u 

<C'^|(/>|o^ _1 2 |V(/>| 2 |(/>div(9^u| 2 + |^>|o^ _1 2 | V 2 (/>| 2 |(/>div9^u| 2 
+ |0|o^ _1 3 |V^| 6 |V</)|3|^div5^| 2 ^ 

<c(\ct>\^ 2 \\ V^|| 2 + |0|^^” 3 |V0| 6 |V0| 3 ) |0V(^«)| 2 
<^ 0 ^ |^V(^«)| 2 

47 - 2 , 5-2 _ 

<Cc 0 5-1 + -iV^ + ^v^n, 

which, combining with (13.261) and 5 > 1, implies that 

h < ^|>/*^V(fl£«)|! + C|u| 2 + Car, for |C| < 2. 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


(3.28) 
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For the term 1 . 4, it is not difficult to show that 

h = J V0 2 • diQ{y) ■ diu 

<C , | 0 |oo|V^| oo ||V«|| 2 |^«| 2 < C4\diu\ 2 , for ICI < 2. 
For the term 1 $, we have 

h = ~ j (d$(v ■ Vu) - v • V(d£u)) • diu 
<|^(u • Vu) - v ■ V(^«)| 2 |^«| 2 
^(IVvIooIIVuIIi + |^| 6 |V«| 3 )|^«| 2 
<Cc 3 \\u\\l + Cc 3 \Vu\ 3 \diu \2 < C'c 3 ||«|||, for 1 < |C| < 2. 

For the term Iq, when |d = 1 we have, 

h = ~ J (di{(ct) 2 + rj 2 )Lu) - (0 2 + ri 2 )Ldiu) ■ d$u 

<C\di{((t) 2 + if)Lu) - (^) 2 + r] 2 )Ldiu\ 2 \diu\ 2 
<C\V((j) 2 + 77 2 )| 00 |iit| 2 |^it| 2 
<C|0|oo|V0|oo|ubiMl>a < C'cgllV^llf. 

When |d = 2, we have 

i e = -J{o^ + e)Lu)-^ + ^)LSiu)-aiu 
<C J ^\\7 2 4>\\(j)Lu\ + | V(/)| 2 |Lu| + \4>VLu\\X7\diu\ 
<c(|V 2 0| 3 |0V 2 u| 6 |V 2 u| 2 + |V<CMd 2 + iv^ioo^v^biv 2 !^) 

<c||v<^|| 2 K 2 + ^l^v 3 dl 

<Ccl\u \ 2 D 2 + + r/ 2 V 3 u||, 

where we used the fact that 

|0V 2 tx| 6 < C\cpV 2 u\ D i < C(|Vdoo|V 2 u| 2 + |</»V 3 «| 2 ). 

Combining (l3.31H - d3.32p . we then have 

h<^\V^T^V 3 u\ 2 2 + Cc 2 0 \\Vu\\j, for 1 < |C| < 2. 
Term Ij can be controlled as follows, 

h = j (a|(V0 2 • Q(v)) - V«/> 2 • diQ(v)) • diu 


(3.29) 


(3.30) 


(3.31) 


(3.32) 


(3.33) 


(3.34) 


<\di(V(j) 2 Q(v)) - V</) 2 a|Q(u)| 2 |^u| 2 . 


(3.35) 
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When ICI = 1, 

l^x(V0 2 • Q(v)) — V0 2 • dfcQ{v )\2 

+ l^loolVVblVvU) < Cel 

when |C| = 2, 

\di{V4> 2 -Q(v))-Vct> 2 -diQ{v)\ 2 

<c(|V<C|V 2 u | 2 + |^|oo|V 2 0 | 6 |V 2 ^| 3 

+ |0|oo|V 3 0| 2 |Vv|oo + |V0| O o|V 2 0| 2 |V77| oo ) 

<Cc|, 

which, combining with (|3.35lh implies that 

h <Ccl\d$u\ 2 , for 1 < |C| < 2. 

Therefore, from (I3.23D - (I3.24D . (I3.28I1 - (I3.30D . (|3.34l) and d3.38ll . we deduce that 

^IMli + I + ?? 2 V 3 u|l < Ccj\\u\\l + Ccf. 

Now Gronwall’s inequality implies that 

IKi)||i+ [ W<t> 2 + rj 2 V 3 u\lds < (Ccgt + ||it 0 ||i) exp(C'cgi) < Ccg, 

Jo 

for 0 < t < T 2 = min (T*, (1 + c^)~ 2K ). 


(3.36) 


(3.37) 


(3.38) 

(3.39) 


(3.40) 


Step 2. We estimate \dx.ut \2 when |C| < 1. From the momentum equations d3.7l) 9. we have 

A'y 27-2 


I fob = 


v • Vu + 


7 - 1 


VC>— + (^ + n z )Lu - V(/) z • Q(u) 


27-2 t 

<c (|u| 6 |Vu | 3 + IC’IcCt 1 |V0| 2 + |</> 2 + 7 2 |oo|«Ie) 2 + |C>|oo|vC)|oo|Vu | 2 


<Ccf. 

Similarly, for |it* | ^»i, we have 

Mr>* = 


Ary 2'y —2 

v • Vu -|- \7<j) 1 + (0 2 + if)Lu — V</> 2 • Q(v) 

7- 1 


D 1 

<C( ||u|| 2 ||Vu||i + |0|^ _1 |V 2 C»| 2 + 


+ |\/V> 2 + ?? 2 |oo|\/^ 2 + 7 2 V 3 u| 2 + |C>|oo|VC>|ooklD 2 


+ MI + l|vc>|||)||Vu|h 


<c(cf+ c 0 |a/C> 2 + 7 2 V 3 u| 2 ), 


which implies that 


/ Kl^ids < C [ (cl K +cl\V¥T^V\\l)ds < Ccg, 

Jo Jo 


(3.41) 


(3.42) 
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for 0 < t < T 2 = min(T*, (1 + C3) 2K ). 


□ 


Finally, we give the higher order estimates for velocity u such as \u\l°°([o,t];D 3 ) and 
|</>V 4 'u|£ 2 ([ 0 j t] ; £ 2 ) through an accurate analysis on the artificial viscosity (f 2 + i] 2 )Lu. 

Lemma 3.4. Let ( <j>,u)(x,t ) be the unique strong solution to \3. 7| ) on M 3 x [0,T]. Then 

\u(t)\ 2 D 3 + \u t {t)\ 2 D i + j (\cj)V 4 u{s)\l + \ut(s)\ 2 D 2 ^ ds < Cc 2 2 K 

for 0 < t < T 3 = min(T*, (1 + C 3 )“ 2A ~ 4 ). 

Proof. The proof is divided into two steps. 

Step 1. The estimate of \diu \2 for |£| = 3. Multiplying (13.221) by diu on both sides and 
integrating over K 3 by parts, we get 


1 d 


- — \d^u\l + a\\Z(j) 2 + r] 2 'V(d^.u)\ 2 + (a + f3)\\/ <f> 2 + r/ 2 div9|u| 


2 dt 


- j{v -V(d^u)) ■ d^u - 5 —^~ J (V(0 2 + r, 2 ) ■ Q(df>u)) ■ d<u 

- J (^((^> 2 + V 2 )Lu) - (<f 2 + r) 2 )Ldiuj ■ diu 

Ea 


+ 

14 


(3.43) 


Now we estimate Ij (i = 8 , • • • , 14) term by term. Similarly to the derivation of (13.241) . 
the first two terms can be treated as follows, 


h 


/ 


(v ■ V(d£u)) ■ d^u < Ccs\d^.u 


2 
2 ’ 


5-1 


(V (cj) 2 + v 2 ) ■ Q(diu)) • diu < -|V^+^V(^«)|1 + Cc 2 1 


lull 


(3.44) 


h = ~ 


5 
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For the term / 10 , when |£| = 3, one has 


ho = - -^-r [ ■ diu = -^-L [ dU^divdiu 



u 


7-1 J x 7-1 

<cj (|0|^- 4 |V0| 3 + |^|^r- 3 |V0||V 2 0| H- |^|^- 2 |V 3 0|)|0di V ^ , 

<c(^VMVct>\U^~^V{diu )\ 2 + |V 2 0| 2 |V0 |oo|(/)|J^” 3 |(/>V(^m )| 2 

2 ) 

4 7-4 ^ 


+ |v 3 «/»| 2 |«/»u 2 


5 1 2 |0V(a£u)|: 


(|V0|H 8 + |V^| 2 |V^hF 6 + |V 3 0|||0|^ 


47 — 4 


<c 


a 

20 
47 - 26-2 


+ -|</>V(d ^)| 2 


<Cc 0 ^ 


+ —IV^ 2 + 7 2 V(9^)||. 


(3.45) 


For the term hi, from integration by parts, we have 


In = j V 0 2 • d$.Q(v) ■ d^u 

<c j (|v 2 0||v 3 ^||0^^| + |v0| 2 |v 3 ^||oS^| + |v0||v 3 ^||0va£^|) 

<c(|V 2 0| 3 |V 3 i ;| 2 |^| 6 + |V</)|^|V 3 i ;| 2 |^| 2 + |V</>UV 3 i;| 2 |0V(c^)| 2 ) ( 3 ‘ 46 ) 

<c(|V 2 0| 3 |V 3 i;| 2 (|^V(9 ^)| 2 + |V^|oo|^«| 2 ) + c||^u | 2 + c 2 |<)>V(^)| 2 ) 

<c(4\d£u\ 2 + 4) + ^\V^ 2 + ifv(d^u)\l, 


where we used the fact that 


|0V 3 xi | 6 < C\c^\i\ D i < C(|V0UIV 3 u | 2 + \^u\ 2 ) . (3.47) 

For the term I 12 , letting r = b = 2 , a = + 00 , / = v, g = Vtt in (12.31) of Lemma [2731 
from (13.401) we obtain 

A 2 =-/(^.V«)-,.V(^)).^ 

<|^(u • Vu) - v • V (d^u)\ 2 \d^u\ 2 

<C (|Vu|oo|V 3 ii | 2 + |V 3 u| 2 |Vu|oo) \diu\ 2 (-3 48 ) 

<C(c 3 |V 3 u | 2 + c 3 ||Vu|| 2 |V 3 k| 2 ) 

<C (c 3 |V 3 u | 2 + c 3 (c 0 + |V 3 u| 2 )|V 3 n| 2 ) 

<C(c 2 |V 3 u | 2 + c 2 ). 
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For the term 1 13 , from |£| = 3 and (13.471) we have 


Ii3 = ~ J (3 + rf)Lu) - (0 2 + rf)Ldiu) • d c x u 

<c j (iWlMI^I + |v*nv>*|M|fi£«| + |v*|»|vl«||*»|) 

+ c/(|vVll^|| 8 c„| + |,v^||v^||^|) 

<c(|V 3 0| 2 |V 2 u| 3 |</>V 3 u| 6 + |V0 |oo|V 2 0| 3 |V 2 w| 6 |V 3 ^|2 + |V0|L|V\|| 
+ |V 2 ^| 3 |V 3 7X| 2 |0V 3 7X| 6 + |V0|oo|0V 4 u| 2 |V 3 u| 2 ) 
<C^c 0 |V 2 'u|||V 3 tx|| (|0V 4 «| 2 + |v</>|oo|v 3 u| 2 ) 

+ CoMl>3(|</>V 4 u| 2 + |V</>UV 3 u| 2 ) + Colrtl^a) + ^\(j)V 4 u\l 

<c (4Mb + <% K+A ) + ^\V<t> 2 + v 2 v 4 u\l 


For the term 1 14 , we notice that 


d£(V</> 2 -Q(u))-V</> 2 -^Q(u) 

= E l m (cujkVd^ 2 • + c^va^y • Q(u)) + v<9^ 2 • q(u), 

l<i,j,k<3 


where £ = C 1 + C 2 + C 3 for three multi-indexes ( l £ R 3 (i = 1,2,3) satisfying |£ l | = 1 ; 
Cujk and C‘ 2 ijk are all constants. The summation is through all permutation on i,j,k. 
hjk = 1 if z, j and k are different from each other, and otherwise Ujk = 0 since there is no 
duplicated differentiation with respect to such decomposition on (. Then, we deduce that 


h* = J (d£(V0 2 • Q(v)) - Vc/) 2 ■ d$.Q(v)) ■ d Q x u 

= J ( E 'ijkCiijkVtfb 4> 2 ■ dx Hk Q{v)) ■ diu 

i,j,k 

+ j ( E l ijkC 2 ijkVd£+< k (j) 2 • d£ Q(v)) ■ diu + 

i,j,k 

= :/l41 + II 42 + 4 143 . 


(3.50) 

(V«9 ^ 2 • Q(v)) ■ diu 
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It follows that 

I\A\ = f (^kjkCiijkS?(t> 2 d c x J+C ’ k Qiy)^ ■ diu 
i,j,k 

<C (|V0|^|V 3 n| 2 |V 3 i;| 2 + |V 3 t>| 2 |0V 3 U | 6 |V 2 </>| 3 ) 

<C (c||V 3 u| 2 + c|(|\/ (j> 2 + r/ 2 V 4 u| 2 + |V</>|oo|V 3 n| 2 )) 
<C (c||tt|£)3 + 4) + — | \/4> 2 + 7? 2 V 4 k|2, 

I\42 = /(E hjkC 2 ijkVd %' +?fc ^> 2 • d^Q(v)j ■ cf x u 
i,j,k 

<C (|V 3 ^| 2 |V 2 ^| 3 |^V 3 ^| 6 + |V0| oo |V 2 0| 3 |V 2 i;| 6 |V 3 ii| 2 ) 
<C (c||V 3 a| 2 + c 2 (|i/ (j) 2 + t/ 2 V 4 u| 2 + |V^|oo|V 3 u| 2 )) 

<C(c 3 H D 3 +c|) + ^|^ 2 + ?? 2 V 4 u| 2 . 

For the term /i 4 3 , from integration by parts, we have 


(3.51) 


^143 



■diu 


= - f Y1 i 9 * c V< ^ 2 ■ d xQ( y ) • d x u ) - j ( d ? Cv ^ 2 • Q ( v ) • ^ +c ‘«) ^ 3 ' 52 ^ 

i=l i =1 

= 'Il431 + ^1432- 


For simplicity, we only consider the case that i = 1, the rest terms can be estimated 
similarly. When i = 1, the corresponding term in /i 4 3 i is 

J i 43 i = - J d x ^V ^ 2 • dfQ(v) • 

<C (|V 3 ^>| 2 |V 2 i;| 3 |(/>V 3 «|6 + |V^»| oo |V 2 0| 6 |V 2 i;| 3 |V 3 n| 2 ) (3 g3) 

<C (c|M D 3 + c|(|0V 4 tt| 2 + |V</>V 3 u| 2 )) 

<C {4\u\ D 3+4) + ^W4> 2 +rt 2 \7 A u\l, 

and the corresponding term in /i 4 32 is 



=:Ia + Ib + Ic + Id- 


(3.54) 
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It is not hard to show that 


/, = - 2 j (af^v*). gw- 


<c , |v 3 <^| 2 |Vv|o 0 |<^v 4 «| 2 < c 4 + + 

For the term Ib, from integration by parts, we deduce that 


A u\l 


Ib = — 2 
= 2 


J (px • Q(u) • 9S +Cl u 

/ < 9 ^ (dfVc!)df(i>)^ -Q(v) ■ diu 
<C J (jV 3 u||Vu|(|V 2 0| 2 + |V^||V 3 0|) + |V 3 u||V 2 u||V 
<C'(|V 2 ^| 3 |V 2 (/ ) | 6 |V 3 u| 2 |Vu| 0O + |V 3 u| 2 |Vu| oo |V 3 0| 2 |V ( / ) | 
+ |V 3 u| 2 |V 2 u|3|V 2 «/»| 6 |V«/»| 0O ) 


<Ccl|V 3 u| 


2- 


Similarly to Is, we have 

-fc + Id < C'c||V 3 ri| 2 , 

which, combined with fl3.50D - (13.57jl . implies that 


hi < —1\/^ 2 + ?y 2 V 4 u|| + Cc3|u| D 3 + (7 c 4 . 


r,2K + 4 

^3 


So, from (|3.43j) - (j3.49|) and (|3.58p , we arrive at 

^\ u \d3 + f (<t> 2 + ?? 2 )|V 4 u| 2 < C4\u\ 2 d3 + (7 c 2 

Then from Gronwall’s inequality, we have 

\u(t )\ 2 D 3 + [ |\/</> 2 + t/ 2 V 4 m|| ds < (7(|u 0 |£)3 + c 2A+4 t) exp((7cff) < Ccq, 
Vo 

for 0 < t < T 3 = min{T*, (1 + i) _2A_4 }. 

Step 2. The estimate for \diut \2 when 1 < |£| < 2. First, from (13.421) we have 


(3.55) 


(3.56) 


(3.57) 


(3.58) 


(3.59) 


(3.60) 


Mdi <C (cf + c 0 |y/(/> 2 + r/ 2 V 3 u| 2 ) <C(cf+ Cq) < Cc 


R 

2 ■ 


(3.61) 
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Similarly, from the momentum equations we also have 


\ut\o 2 — 


Ary 2 ^ — 2 

v ■ Vrt H-V0 S - 1 + (cf) 2 + r] 2 )Lu — V0 2 ■ Q(v) 

7-1 


<C( |M| 3 ||Vu|| 2 + Ml-* V 3 ^| 2 + 


D 2 

2 , 


2-/-2 


|V0| 3 |V 2 0| 6 

+ + Iv^+^loolv^+^v^b 

+ (|0|oo + ||V(/)||2) 2 (||rt||3 + IMb)) 

<C (cf + colv^T^V^b) , 


which implies that 

rT 3 


J \u t \Us< C f (c 2K + c 2 \V^+V 2 ^u\t) ds<Cci 


Combining the estimates obtained in Lemmas 13.2113.41 we have 

i + 4 > + |0(f)& + ||0(t) — 0H3 <Ccq, 

< Ccf, \<fh{t)\ 2 D i < Cel, \M0 \d 2 <Ccl 

IK*)Hi < Cel, \u(t)\jj 2 + \ut(t)\l + J (|0V 3 u|l + \u t \ 2 D Jjds <Cc\ K , 

\ u (^)\‘d3 + \ u t(t)\i)i + J (| 0 V 4 it |2 + |rit||) 2 ^ds <Cc 2K , 

for 0 < f < min(T*, (1 + C 3 ) _2 ^~ 4 ). Therefore, if we define the constants c, (i = 
and T* by 


ci = C 2 c 0 , 

c 4 = C2c| = 


c 2 = C^cf = , c 3 = ch« = C ££i ^ ±i cf, 

C x2+K+ §c^ 2 , and T* = min(T, (1 + c 3 )~ 2K ~ 4 ), 


then we conclude that 


l + 0 2 + sup (| 0 (f)|^ + || 0 (f) - 0 ||! + |K*)||i) <cf, 
0 <t<T* 

sup (|0t(f)|i + |w(t )||»2 + Mf)||) + [ (\<l>V 3 u\l + \ut\ 2 Dl )dt <c|, 

0 <£<T* Jo K ' 

rT* 

sup (|0t(f)|^i + |tt(f )||,3 + Mf)||,i) + / (| 0 V 4 u|| + K|^ 2 )df <c|, 

<t<T* J 0 v 7 

sup | 0 t(f )||,2 <cf 

0 <t<T* 


(3.62) 

□ 

(3.63) 

1 , 2 , 3 , 4 ) 

(3.64) 


ess 


(3.65) 
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3.4. Passing to the limit 77 —>• 0. With the help of the ^-independent estimates estab¬ 
lished in ()3.65p . we now establish the local existence result for the following linearized 
problem without artificial viscosity under the assumption (f >0 > 0 , 

f £ -y 

4>t + v ■ V(f -I-— 'Ipdivv = 0, 

u t + v-\7u -1- 'l-V(j)ir=T +(j) 2 Lu = V0 2 • Q(v), 

< 7—1 (3.66) 

0, u)\ t =o = (4 >q(x),u q (x)), x € M 3 , 

((f), u) —>• (<f, 0 ), as \x\ —v + 00 , t > 0 . 

Lemma 3.5. Assume that the initial data satisfy 113. 5)) . Then there exists a unique strong 
solution ( 1 <f>,u)(x,t ) to 13. <?£>!) such that 

cf)-f€ C([0, T*]-,H 3 ), eft € C([0, T*\,H 2 ), 

ueC([0,T*y,H s ')nL°°([0,T*];H 3 ), (fV 4 u € L 2 ([0, T*]; L 2 ), (3.67) 

€ C([ 0 , T*]; 77 1 ) nL 2 ([ 0 ,T*]; Z7 2 ) 

for any constant s' € [2,3). Moreover, (<f,u) also satisfies the a priori estimates in L3. 65\) . 
Proof. We shall prove the existence, uniqueness and time continuity in three steps. 


Step 1. Existence. From Lemma [3.11 for every 7 > 0, there exists a unique strong solution 
((fp, u v )(x, t) to the linearized problem (13.71) satisfying the estimates in (13.651) . which are 
independent of the artificial viscosity coefficient 7 . 

By virtue or the uniform estimates in (13.651) independent of 77 and the compactness in 
Lemma 12.21 (see ESI), we know that for any R > 0, there exists a subsequence of solutions 
(still denoted by) ((f> ri ,u ri ), which converges to a limit ((f), u) in the following strong sense: 

-> ((f),u) in C([0,T*]-H 2 (B r )), as 77 -> 0. (3.68) 


Again, using the uniform estimates in (13.651) independent of 77 , we also know that there 
exists a subsequence (of subsequence chosen above) of solutions (still denoted by) ((fP, u v ), 
which converges to ((f), u) in the following weak or weak* sense: 


(0 V 


((f), u) 

weakly* 

in 

L°°([0,T*]-H 3 (R 3 )), 

(ft 

weakly* 

in 

L°°([0,T*]-,H 2 (R 3 )), 

Ut 

weakly* 

in 

L°°([ 0 , T*]; Lf 1 ^ 3 )), 

Ut 

weakly 

in 

L 2 ([0,r*];L> 2 (M 3 )). 


(3.69) 


Combining the strong convergence in (13.681) and the weak convergence in (13.691) . we 
easily obtain that ((f), u) also satisfies the local estimates in (13.651) and 


<^VV -»• </>V 4 u weakly in L 2 (R 3 x [0,T*]). 


(3.70) 


Now we are going to show that ((f), u) is a weak solution in the sense of distribution to 
the linearized problem (13.661) . Multiplying (13.71) 9 by test function w(t,x ) = (w 1 ,w 2 ,vj 3 ) € 
C'^°(M 3 x [0, T*)) on both sides, and integrating over K 3 x [0, T*], we have 
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/ 0 


^ ■ wt — (v ■ V)u v ■ w 4—divu;) dxds 


= - / uq ■ w(0,x) + 


I J (^((fa) 2 + t] 2 )Lu v ■ w — S/fa 11 ) 2 ■ Q(v) ■ w^jdxds. 


(3.71) 


Combining the strong convergence in (13.6811 and the weak convergences in (|3.69l) - (l3.70p . 
and letting r/ —> 0 in (13.711) , we have 



(3.72) 


A'y 2 7 -2 \ 

u ■ wt — (v • V)u • w H- 4> *- 1 divu> dxds 

7 — 1 / 

= —J uq-w(0,x) + J J i^) 2 Lu-w— Vcf) 2 ■ Q{v) ■ vo^jdxds. 

So it is obvious that (0, u) is a weak solution in the sense of distribution to the linearized 
problem (13.661) . satisfying the following regularities 

0 - 0 € L°°([0, T*]-H 3 ), fa € L°°([0, T*]-,H 2 ), 

ue L°°{[0,T*]-H 3 ), faS/^u € L 2 ([0, T*]; L 2 ), (3.73) 

u t € L°°([0, T*]; H 1 ) n L 2 ([0, T*]; D 2 ), 

where we used the lower semi-continuity of various norms in the weak or weak* conver¬ 
gence in (I3.69I) - (I3.70I) . Therefore, this weak solutions (0, u) of (13.661) is actually a strong 
one. 


Step 2. Uniqueness. Let (fa,ui) and (0 2 ,u 2 ) be two solutions obtained in the above step. 
For ip = fa — fa, we have from (13.661) i that 


<Pt + v ■ V<p = 0, (3.74) 

which immediately implies that ip = 0 in R 3 with zero initial data. 

For u = ui — U 2 , from (13.661) 9. using the fact fa = fa, it is clear that 

Tit + v ■ Vu — (film = 0. (3.75) 

Multiplying (13.751) by u on both sides, and integrating over R 3 , we have 

4l\u\1 + l^iVu|| <C\Vv\oo\u\l + C , |u| 2 |V0i| oo |0iVu| 2 
at (3.76) 

<-^\faVu\l + CdVuU + \Vfa\D\u\l 

Now, the Gronwall’s inequality, along with zero initial data of u implies that u = 0 in R 3 . 
This completes the proof of uniqueness. 

Step 3. Time continuity. For 0, the regularities in (13.731) and the classical Sobolev imbed¬ 
ding theorem infer that 


0 - 0 € C([0, T*];H 2 ) n C([0, T*]; weak - H 3 ). 


(3.77) 
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Using the same arguments as in the proof of Lemma 13.21 we have 




<(H0o - 0\\l + c y o (IIW'llllMl! + \<pv 4 v\l)ds) exp [C (||u(s)|| 3 + l)dsj, 

which implies that 


lim sup 
t-to 


b < 


(3.78) 


(3.79) 


From Lemma 12.71 and (13.771) . we know that </> is right continuous at t = 0 in H 3 space. 
The time reversibility of the equation (13.661) i yields 

</>-^eC([0,T*];i7 3 ). 

For (j>t, we note that 


cf) t = —v ■ Vcj) ——— ipdivv. 


(3.80) 

(3.81) 

(3.82) 

(3.83) 


On the other hand, since 

V’Vu € L 2 ([0,T*];i7 3 ), (^Vv) t € L 2 ([0, T*]; H 1 ), 
using Sobolev imbedding theorem, we have 

4>VveC([0,T*]-,H 2 ), 

which implies that 

4>t £ C([0, T*]-, H 2 ). 

For velocity u, from the regularities shown in (13.731) and Sobolev imbedding theorem, 
we know that 

u € C([0, T*]-H 2 ) 0 C([0, T*]; weak - H 3 ). 

Then from Lemma 12.51 for any s' € [2, 3), we have 


(3.84) 


IMIs'< C3IMI0 3 \\ u \\i • 

This, together with the upper bound estimates shown in (13.651) and the time continuity 
(USD, yields 


ueC([0,T*];H s ). 
Finally, we consider ut ■ Noting that 


2 Ary 2r—5 _1 

ut = —v ■ Vu — -— -<t> S ~ 1 Vcj) + (j?Lu + Vc /) 2 ■ Q(v), 
o — l 


(3.85) 

(3.86) 


where 


€ L 2 ([0,T*];H 2 ), 


Q ( v ) = — [a(Vv + (Vv) T ) + 
we then have from (13.731) that 

4> 2 Lu € L 2 ([0, T*]; H 2 ), (< j> 2 Lu) t € L 2 ([0, T*]; L 2 ), 

which means that 

<; b 2 Lu G C([0, T*]; H 1 ). 


(3.87) 

(3.88) 
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Combining (13.91) . (13.801) . (I3.85[) and (13.88[) . we deduce that 

Ut G C([0, T*];!! 1 ). 

We thus complete the proof of this lemma. 

□ 

3.5. Proof of Theorem 13.11 Now we turn to the nonlinear problem (I3.1I) - (I3.4I) . Our 
proof is based on the classical iteration scheme and the existence results for the linearized 
problem obtained in subsection 3.4. We first define constants Co and ci, C2, C3, C4 as in 
subsection 3.3. Assume that 


2 + (t> + |0o 100 + ||(0o - 0) ^0) 11 3 < c 0 . 

Let (0°,Ti°) be the solution to the system 

Y t + uq ■ VT = 0 in (0, +00) x M 3 , 

Zt — Y 2 AZ = 0 in (0, +00) x R 3 , 

(Y,Z)\ t=0 = ( 0 o, no) in 

(T, Z) -A- (0,0) as |x| —> +00, t > 0, 

with the regularities 

0°-0GC([O,T*];fL 3 ), 0°V 4 u° £ f 2 ([0,f ];I 2 ), 

u° G C{[0,T*]-,H s ')nL°°([0,T*]-H 3 ) for any s'€ [2,3). 

Due to the regularity of (0°,ti°)(x), there is a positive time T** G (0, T*] such that 
1 + 0 2 + _ sup (|0°(t)|L + 110° (t) - 0lli + ll' u °(^)lli) <ci, 

rT* 


(3.89) 


(3.90) 


0 <t<T 

SUp (|0?(t)[| + W°(t)\ 2 d2 + |li?(t)H) + 


0 <t<T 


ess 


V 3 u°|2 + \Uf ||)i^dt <4, 

sup m 2 D1 +\u% 3 +\^i 2 D1 )(t)+ r (i0°wii+k% 2 )dt< c |, 

0 <t<T** Jo ' 7 

SUP 10? 11,2 <4 
0 <t<T** 

We now give the proof of Theorem 13.11 

Proof. We prove the existence, uniqueness and time continuity in three steps. 


(3.91) 


Step 1: Existence. Letting (if,v) = (0°,u°), we first define (0 1 ,ti 1 ) as a strong solution 
to problem (13.661) . Then we construct approximate solutions (<f> k+1 ,u k+1 ) inductively, as 
follows: assuming that (<j) k ,u k ) was defined for k > 1, let (0 fc+1 ,u fc+1 ) be the unique 
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solution to problem (13.66[) with (^,v) replaced by ( 4> k ,u k ), i.e., (4> k+1 ,u k+1 ) is the unique 
solution of the following problem 

$+1 + u k . v ^+1 + iZ-I/dm/' = 0j 

u k t +l + u k ■ Vu k+1 + |^-$ fc+ 1 V </> fe+1 + (4b k+l ) 2 Lu k+l = V(c^ +1 ) 2 • Q(u k ), 
((f> k+ 1 ,u k+ 1 ) |t = 0 = (</>o,«o), x € M 3 , 

(</> fc+1 , rt fc+1 ) —>• (0, 0), as |cc| —>• +oo, t > 0, 

where <f> fc+1 = (. cf > k+1 ) «-i . 

From the estimates shown in Section 3.4, we know that the sequence ( 4> k ,u k ) satisfies 
the uniform a priori estimates in (13.651) for 0 < t < T**. 

Now we prove the convergence of the whole sequence ( 4> k ,u k ) of approximate solutions 
to a limit (</), u) in some strong sense. Let 

0 k+1 = ( j) k + 1 -0 k , u k+l = u k+l - u k , 

$ fc+1 = $ fc+1 - <]? fc = (0 fc+1 ) 21 s-i 1 _ (0 fc ) 2 T-i~ 1 . 

Then from (13.921) we have 

^' +1 + u k ■ Vf +1 + u k • V0 fc + + (j) k divu k ) = 0, 

u k+1 + u k ■ Vu k+1 + u k ■ Vu k + |^l($ fc+1 V^ fc+1 + ¥ fc+1 V^ fc ) 

< d - 1 (3.93) 

= -(4> k+1 fLu k+l -4> k+1 ((j) k+1 + cj) k )Lu k 

, +V (4> k+1 (<j> k+1 + 0 fc )) • Q(u fc ) + V(<^) 2 • (Q(u fc ) - Qfa*" 1 )). 

First, multiplying (13.931) i by 2cj) k+1 and integrating over R 3 , we have 

A|/ +1 |2 = _ 2 j [u k ■ V/ +1 + u k ■ V(t> k + (^ fe divu fc_1 + <^dml fc ))^ fc+1 
<C'(^|Vn fc | 00 |^ fc+1 || + |^ +1 |2|^|2|V/|oo 

+ |/| 2 |V« fc - 1 |oo |/ +1 |2 + |^ fc divu fc | 2 |</> fc+ 1 | 2 ), 

which yields that ( for 0 < v < ^ is a constant) 

’4|/ +1 («)|| < ^(t)i/ +1 (t)i| + -(i^wil + i/(t)H + i/div«‘(t)H), 

< 

A k {t) = c( |V« fc |oo + —|Vu fc_1 |^ 0 + -|V0 fc |^ + -). 

\ V V V / 


(3.94) 
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From (13.651) . we know 


A k (s)ds < C v t, for te[0,T**}, 


where C v is a positive constant depending on u and constant C. 
Second, we multiply (13.931) 9 by 2 u k+1 and integrate over M 3 to find 

d 'u k+1 \l + 2a\4> k+1 Vu k+1 \l + 2 (a + f3)\c/) k+1 dwu k+1 \l 


dt 

= - 2 


u k • X7u k+1 + u k • W 


k\ „-h+1 


- 2 


2^7 /^fc+iv 7 j; fc + 1 , Tf7 fc+1 v7A^A tA+i 


<5-1 


($ fc+i V</> " + $ " V0 fc ) • if 


-4 a J (t> k+l V(j) k+l ■ \7u k+1 • if +1 - 4(a + P) J f fc+1 V</> fc+1 • U fc+1 divlf +1 


- 2 


-fc+i, fc+1 


+ 0 fc )Lu* • u* +i + 2 J V(/' +1 (0 fc+1 + </> fc )) • Q(u fc ) • u 

9 


■ fc +1 


+ 2 / V(f fc ) 2 • (Q(u fc ) - Q^- 1 )) ■ u k+1 =: Y, J i ■ 


2=1 


(3.95) 

We now estimate Jj (i = 1, ■ ■ ■ , 9) term by term. For the term J), we see from integra¬ 
tion by parts that 


Ji = -2 J u k ■ Vu k+1 -u k+1 < C\S7u k \ 00 \u k+l \ 2 2 . 
For J2, it is easy to show that 

J 2 = - 2 [ u k ■ \7u k ■ u k+1 


(3.96) 


(3.97) 


<C|Vu fc | 0O |F fc | 2 |F fc+1 |2 < -\Vu k \lJ[u k+l \ 2 + v\u k \ 2 2 . 
Applying integration by parts for J3, we have 
J3 = — 2 


Ml^yf + 1 -u k+1 

5 — 1 


4A7 

= J- 

<C 


<C 


l J ^+^ k+1 dwu k+1 + t +1 ™ k+1 -^ k+1 




• u 


l ) 


fe+1 | 2 - f |/ +1 | 2 |^ +1 divu fc+1 | 2 + |/ +1 | 2 |V$ fe+1 |oo|lZ fe+1 | 2 

k+1 \M^\t +1 \l + \t +1 \l + |v^+ 1 |i 0 |^-+ 1 |l) + ^ fc +W +1 | 2 . 


(3.98) 


J4 is estimated directly below, 


J 4 = - 2 


Tfi fc + 1 V7^fc . 777^ + 1 


[Ml 

J 5-1 


$ " V(/) K • u K 


<C|4> fc+1 | 2 |V^ 


A; I I7J7A7+I1 
00 U 2 


(3.99) 


<C 


fc+ii 


\<t> k 


tM-fc+1,2 


5-1 
loo 




k |2 I—fc-(-l |2 


m 


12 > 
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where we used the fact that 

l ^' +1 |2 = |(^ +1 )^ - (<^)“| 2 < C\\cj> k+1 \ + 
Similarly, we are able to treat terms J 5 -J 7 in the following way, 

J 5 = - 4a I (j) k+l \7(t> k+l ■ Vu k+l ■ u k+l 

<C\V(t> k+l \^\(t> k+l Vu k+ %\u k+1 \ 2 

<C|V 0 fc+ 1 |Ll^ +1 || + |^ fe +W +1 ||, 

J 6 =- 4 (a + P) J cf +l V(j ) k+1 ■ u k+ 1 dwu k+1 

<C|V0 fc+1 |oo|0 fc+1 Vn fc+1 |2|u fc+1 |2 

<C|V ^ +1 |^ fe+1 | 2 + W +1 \l 

J 7 =2 J ^ k+ \(j) k+1 + (f> k )Lu k • a k+1 


^1 (j) k+1 

00 1 ^ 


2- 


<c(|^ +1 |2|</» fc+1 R fc+1 | 6 |Lu fc | 3 + |/ +1 | 2 |^Lu fc |oo|^ +1 | 2 ) 

<c(|^ +1 | 2 (|^ + W +1 | 2 + |V/ +1 | 00 |R fe+1 | 2 )|LR fe | 3 

+ |/ +1 |2(||V0 fc ||2||^||3 + I/'W|2)|R fc+1 | 2 ) 

<c(\Lu k \l\f +l \l + \f +l \l + iV^^lL^lllu^H 

+ (||V0 fc [| 2 || W||! + |</> fe vV| 2 )V +1 |^) + ^ k+ 1 Vu k+1 \l 


(3.100) 


where we used the fact (see Lemma l2.1D that 

l/Wloo <C| 0 fc VV|||V(</. fc VV)|| 

<G|^vV|| 1 |v(^vV )|| 1 
<c ||v(<^ fe vV)||i 

<c , (|V 0 *| oo ||v 2 «*||i + l^lool W | 2 + |<^vV | 2 + |v 2 <^| 6 |vV| 3 ) 
<C , (||V0*|| 2 ||V 2 « fe || 1 + I^VVla). 

(3.101) 


Now we turn to the tricky term J 8 . First we have 

J 8 =2 J V(/ + V +1 + 0 fc )) • Q(u k ) • U fc+1 
= - 2 j J 2^ +1 ^ k+1 + ^ k )di^u k+ ^ - 2 I £ 0 fc+ 1 ( 0 fc+1 + <f)a%djy* +1 ’ j 

i,i i,j 

= : j&l + Js 2 + Js3 + J&4, 


(3.102) 
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where n fc,J represents the j-th component of u k (k > 1 ), 

for k > 1 , j = 1 , 2 ,3, 

and the quantity is given by 


a? = ( a(diU k ’ j + djU k ’ 1 ) + divrt fc <%) for i, j = 1, 2,3, 

where <5^- is the Kronecker symbol satisfying 5ij = 1, i = j, and 5ij = 0, otherwise. 

We are now ready to estimate terms Jsi — Js 4 one by one. First we have 

Jsi = - 2 | ^/ + y +1 d4% fe+1 ^ < C'|vV| 3 |/ +1 | 2 |^ q+ V +1 |6 

<C\V 2 u%\$ k+ \ (|</> fc+1 Vu fc+1 | 2 + |V0 fc+1 |oo|^ fc+1 | 2 ) 

<c (|v</» fc+1 |L|vV|||^ +1 |2 + |n fc+1 || + \V 2 u k \i\t +1 \i) + ^\<t> k+l Vu k+1 \l 

(3.103) 


Similarly, for the terms Is2-h3, using (13.1011) . we get 

J 82 = - 2 j < C'|</> fe VV| 0o |/ +1 | 2 |H fe+1 | 2 

l 1 3 

<C ((llV^llallvVlU + |/W| 2 ) 2 \t +l \l + \u k+1 \l^j , 

J 83 = - 2 J ^0 fe+1 / +1 4diu k+1 j < C'|? +1 | 2 |0 fc+1 V« fc+1 | 2 |Vu fc | oo 

<C\Vu k \l\t +1 \ 2 2 + ^ fe+ W +1 ||. 

For J§ 4 , we have 

^84 = - 2 J 0 k+l 0 k akdiU kJrLj 

= ~ 2 j YA k+1 ^ k - ^ k+1 + 4 +1 )4^ fc+1J 

<C (|Vu fc | 00 |Vu fc+1 | 00 |^ fc+1 |l + |Vir fe | 0O |/ +1 VU fe+1 | 2 |^ +1 | 2 ) 

<c (|V« fc | 00 |V!Z fc+1 | 00 |^ fc+1 |l + |Vn fc |Ll4' +1 ||) + ^|^ +1 VU fe+1 ||. 
Combining ((3.10211 - (13.1051) . we arrive at 

Js <c(|u fc+1 i| + (iv^il + i)i W|§i/ +1 ii + |vu fc iLi4 +1 ll 

+ (||V^|| 2 ||vV||i + 14 W| 2 ) 2 |4 +1 || + |Vu fc | 00 |VH fc+1 | 00 |/ +1 ||) 


(3.104) 


(3.105) 


(3.106) 
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As for the last term Jg, it is easy to show that 

J 9 =2 J V(<^') 2 • ( Q(u k ) - Q(u k ~ 1 )) ■ u k+l 

<C , |V</) fe |oo|0 fc V^| 2 |u fc+1 | 2 < -|V^|^|« fc+1 | 2 + v\ct> k Vu 


(3.107) 


■k 12 
2 > 


In summary, using (13.9711 - (13. 1001) and (13. 106D - (13. 10711 . (13.951) implies that 

d 


— \u k+l \l + a\ct) k+1 Vu k+l \l 
■fe+ii 2 


(3.108) 


<^(i)r +i |2 + I 2 + p(|^Vu fc | 2 + |u fc | 2 ), 

for 0 < v < ^ is a constant. Here B k (t) and B k (t ) are given by 

B k (t) =c(l + |Vu fc |oo + + |Vd> fc+1 |L + + |V0 


1 


fc+l|2 

OO 

Afc |2 
loo 


+ (iiv0 fc ii 2 ||u fe i| 3 + \rv\%y + |v^ +i i^|Lu fc i 2 + -iv0 fc K 

(47 — 45 

1 + |vV|| + ||^ +1 | + |<^||^ + (||v^|| 2 ||u fc || 3 + |0*vV| 2 ) 2 

+ IV^+^IVVI 2 + |Vu fc |L + |Vu fc | 00 |Vu fc+1 | 00 ), 

satisfying the following estimate 

^ (B k (s) + B k {s))ds < C + C v t, t G [ 0 , T**]. 

Jo 


Denote 


r k+1 (t)= sup |/ +1 (s)||+ sup |u fc+1 (s)| 2 - 

se[0,t] sg[0,t] 


From (13.941) and (13.1081) we finally have 

":M k+ \t)\l + \u k+1 {t)\l) + a\ct> k+l Vu k+l \l 


dt 


<E k (|/ +1 (t)|i + |u fc+1 (f)||) + v (\(J k Vu k {t)\l + \cj) k {t)\l + \u k (t)\ 


2 

2 I > 


for some E satisfying 


fE k (s) d. 

Jo 


s < C + C v t. Applying Gronwall’s inequality, we have 


r fc+1 + f a\4>' 

Jo 


k+1 Vu k+1 tids 


<Cv J y\(f) k \7u k \2 + \(f> I 2 + |u fc | 2 Jdsexp (C + C u t) 

<Cvi f \4> k 'Vu k \\ds +1 sup \(j) k (s)\l + t sup |7Z fc (s)|| ) exp(C v t). 
\J 0 se[o,t] se[o,t] J 

We choose v > 0 and T* £ (0, min(l, T**)) small enough such that 

4 Cu < min(a, 1), and exp(C'^T*) < 2. 
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Therefore, we achieve 


J2 (r k+1 (T*) + 


k=1 



a|0 fc+1 Vu fc+1 |ldf) < C < Too. 


(3.109) 


This implies that our approximate solution sequence (c f> k ,u k ) is a Cauchy sequence under 
the topology of L°°([0, T*]; IL 2 (M 3 )), and hence converges to a limit function (0, it) there 
strongly. On the other hand, since the uniform estimates in (13.651) is independent of k. 
we know that there exists a subsequence of solutions (still denoted by) (<j) k ,u k ), which 
converges to a limit (0, u) in the following weak or weak* sense: 


{4> k ,u k ) ->■ ((f), u) 

4>t 0t 
u k Ut 
u k U t 

0 fc VV 0V 4 u 


weakly* in L°°([0, T*]; H 3 {R 3 )), 
weakly* in L°°([0, T*]; H 2 {R 3 )), 
weakly* in L°°([0, T*]; i? 1 (M 3 )), 
weakly in L 2 ([0, T*]; D 2 (M 3 )), 
weakly in f 2 ([0,T*] xl 3 ), 


(3.110) 


which, from the lower semi-continuity of norm for weak or weak* convergence, imply that 
the local estimates in (I3.65P still hold for the limit function (0, it). 

Now, it is easy to show that (0, u) is a weak solution of (I3.1I) - (I3.4D in the sense of 
distribution with the following regularities: 

0-0 G L°°([0,T*];tf 3 ), 0 f € L°°([0, T*\;H 2 ), 

u € L°°([0, T*]; H 5 ), 0V 4 u € L 2 ([0, T*]; L 2 ), (3.111) 

e i°°([o, r*] ; ^r 1 ) n i 2 ([o, .o 2 ). 


The existence of strong solutions is proved. 


Step 2: Uniqueness. Let (0i,«i) and (02,u 2 ) be two strong solutions to Cauchy problem 
(13.11) - (13.41) satisfying the uniform a priori estimates (13.651) . Denote 

^ = 01 - 02, U = Ul - U 2 , 


2 ^ —5 — 1 
f 5-1 


2^y—5 —1 
5-1 

2 


= d>i — $ 2 = 

It follows from (13.11) that ( Jp,u ) satisfies the following system 

<5-1 


(ft + Ul ■ V ip + U ■ V02 + 

Tit + ui • Vu + u ■ Vn 2 + 


2 

2 yly 

<5-1 


(c/?divu 2 + 0idivu) = 0, 

(3>iV^ + <f>V0 2 ) 


= -(0fLrt + <£>(01 + (j)2)LU2) 


+V0 2 (Q(ui) - Q(u 2 )) + V(y?(0l + 0 2 ))Q(lX2), 
with zero initial data. Let 

*(*) = 1^(011 + TOI1- 


(3.112) 
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Using the same arguments as in the derivation of ()3.94|) - (13. 1081) . we have 

^*(t) + C\<kVu(t)\%<F(t)*(t), 


(3.113) 


I F(s)ds < C for 0 < t < T*. 
o 

From Gronwall’s inequality, we conclude that Tp = u = 0. Then the uniqueness is obtained. 


Step 3: Time-continuity. It can be proved by the same arguments as in the proof of 
Lemma 13.51 We omit the details here. 

□ 

3.6. Proof of Theorem 11.11 and Corollary 11.11 Based on Theorem 13.11 we are now 

ready to prove the local existence of regular solution to the original Cauchy problem (11.11) - 

» 


3.6.1. Proof of Theorem 11.11 

Proof. For the initial data (11.161) . we know from Theorem 13.11 that there exists a time 
T* > 0 such that the problem (I3.1I) - (|3.4D admits a unique strong solution (cj),u) satisfying 
the regularities in (|3.6D . which means that 

p^r = (fie C^R 3 x [0,7;]). (3.114) 

Noticing that p = , and > 1 for 1 < <5 < min ^3, TL! j ; it is easy to show that 

p{x,t) G C' 1 (M 3 x [0, T*]). 

Now we verify that (p,u) satisfies the original equations (11.11) . Multiplying both sides 
of (13.11) I by 

do 2 3-5 „ 

v^OM) = € C(R 3 x [ 0 ,T*]), 

we get the continuity equation in ( 11 . 11 ) i . 

Multiplying both sides of (13.11) 9 by 

(f >— 7 = p{x,t) G C 1 (M 3 X [0,T*]), 

we get the momentum equations in ( 11 . 11 ) 9 . Therefore, (p,u) is a solution to (ll.ll) - (U.3D in 
the sense of distribution with the regularities shown in Definition 11.11 
Recalling that p can be represented by the formula 


p(x,t ) = po(U(0] x, t)) exp ^ J divu(s,U(s;x,t))ds^J : 
where U G C 1 ([0,T*] x M 3 x [0, T*]) is the solution to the initial value problem 


£U(s-,x,t) = u(s,U(s;x,t)), 0 < s < T*, 

U(t;x,t) = x, iGl 3 , 0 < t < T*, 


(3.115) 


p(x,t) > 0, V(x,t) G M 3 x [0, T*]. 


it is obvious that 
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In summary, the Cauchy problem (U.1I) - (I1.3D has a unique regular solution ( p,u ). □ 


3.6.2. Proof of Corollary 11.11 . 


Proof. First, from 1 < 6 < - we know that > 3. Since 


4> — (j> £ c([o,t*] ; ft 3 ) n ^([o, n\-,H 2 ), 


and 

p(x,t) = 


we have 

p-p€C([o,uy,H 3 ). 

Second, due to the fact that 

P^r --jrr € C([0,T*];FT 3 ), u € C([0, T*]; H s ') n L°°([0, T*]; FT 3 ), 
€ F 2 ([0, T*]; L 2 ), u t € C([0,T]; 7F 1 ) n L 2 ([0, T]; D 2 ), 


(3.116) 


for any constant s' € [2,3), and the same arguments as in Lemma 13.51 for the time conti¬ 
nuity, we deduce that 

pdivu € L 2 ([0, T*]; F 3 ) n C([0, T*]; 1L 2 ). (3.117) 

At last, with the aid of the continuity equation p t +u-Vp+pdivu = 0 and (13.116I) - (I3.117I) . 
it is clear that 

P - P € C([0, T*];H 3 ) n C\[ 0, T*]; iL 2 ). 

Furthermore, when <5 = 2 and 7 > 3, by the same token, the regularity of p in these 
cases can be achieved. □ 


4. Formation of singularities 

In this section, we consider the formation of singularities of regular solutions obtained 
in Section 3. Two classes of initial data that lead to finite time blow-up will be given. We 
assume that ( p,u)(x,t ) is the regular solution in M 3 x [0 ,T m ) obtained in Theorem 11.11 
with T m the maximal existence time. 

4.1. Blow-up by isolated mass group. The first kind of singularity formation is driven 
by isolated mass group, defined in Definition 11.21 Assume that the initial data (po,uo) 
have an isolated mass group (Aq,Bo), the following definition helps to track the evolution 
of Ao and Bq. 

Definition 4.1 (Particle path and flow map). Let x{t\ xq) be the particle path starting 
from xq at t = 0, i.e., 

^-x(t;x 0 ) = u(x(t;x 0 ),t), x(0]x 0 ) = x 0 . (4.1) 

at 

Let Aft), B(t), B(t ) \ Aft) be the images of A 0 , B 0 , and B 0 \ A 0 , respectively, under the 
flow map of » *-e., 

A ft) = {a:(t;xo)|a;o € A 0 } , 

Bft) = {x(t; x 0 )|x 0 G B 0 } , 

Bft) \ Aft) = {x(t; x 0 )\x 0 € B 0 \ A 0 } . 
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The following lemma confirms the invariance of the volume |^4(t)| for regular solutions. 

Lemma 4.1. Suppose that the initial data (po, uo)(x) have an isolated mass group (Ao, Bq), 
then for the regular solution (p, u)(x, t) on IR 3 x [0, T m ) to the Cauchy problem \1.1\) - H1.3\) . 
we have 

l^(^)l = l-^ol; t £ [0,T m ). 

Proof. Since 

p(x(t; x 0 ),t) = po(x Q ) exp ( J di vu(x(s;x 0 ), s)dsj , 

it is clear that 

p = 0, in Bft) \ Aft). 

From the definition of regular solutions, we have 

ut + u ■ Vu = 0, in B(t)\A(t). 

Therefore, u is invariant along the particle path x(t;xo) with xo € Bo \ Aq. 

For any Xq. Xq E cMq, we define 


^x l (t-,x l 0 ) = u(x l (t,x l 0 ),t), x*(0 ; Xq) = Xq, for * = 1,2. 


Then we have 


d_ 
d t 


(x 1 ^; Xq) - x 2 (t;xl )) = ufx x {t] xj), t) - ufx 2 (t; ®§), t) = 0, 


which implies that 


I^4(^)1 — 1^4-01j t € [0, T m ]. 


(4.2) 

(4.3) 

(4.4) 

□ 


We point out that, although the volume of Aft) is invariant, the vacuum boundary dAft) 
can vary as time evolves. In order to deal with this situation, the following well-known 
Reynolds transport theorem (c.f. HZI) is useful. 

Lemma 4.2. For any G(x,t) € C' 1 (M 3 x M + ), one has 

— ( G(x,t)dx = [ Gt(x,t)dx+ [ G(x,t)(u(x,t) ■ n)dS, 

dt JA(t) JA(t) JdA(t) 


where n is the outward unit normal vector to dAft), and u is the velocity of the fluid. 

In the rest part of this section, we will use the following useful physical quantities on 
the fluids in Aft): 

mft) = / p(x,t )dx (total mass), 

JA(t) 

Mft) = f p(x,t)|x| 2 dx (second moment), 

J A(t) 

Fft) = / p(x, t)u(x, t) ■ xdx (radial component of momentum), 

Jaw i 


eft) = I {-apW\ 2 - -\(x,t)dx (total energy). 

JA(t) v 2 7 - 1 / 
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From the continuity equation, it is clear that the mass is conserved. 

Lemma 4.3. Suppose that the initial data (po, uo)(x) have an isolated mass group (Ao, Bq), 
then for the regular solution (p,u)(x,t) on K 3 x [0, T m ) to the Cauchy problem 
we have 

m(t) = m( 0), for t <E [0,T m ). 

Proof. From (II.lib and Lemma 14.21 direct computation shows 


— m(t ) = f pt dx + f pu ■ n dS 
dt J A(t) JdA(t) 


d A(t) 

which implies that m(t) = m(0 ). 


dA(t) 
—di y(pu) dx = 


Id Apt) 


—pu ■ h d£ = 0, 


□ 


Motivated by m, we define the following functional 

I{t) = M{t) - 2 (t + 1 )F(t) + 2(t + l) 2 e(f) 

= / \x - (t + l)u\ 2 pdx H- ^—(t + l) 5 

JA(t) 7 - 1 


P dx. 


(4.5) 


lA(t) 


We now follow the arguments of m with some proper modifications to prove Theorem 
O One of the key observations in the following proof is that the viscosity tensor T = 0 
in vacuum region due to (11.511 . 

Proof of Theorem 1.2. From system (11.11) . it is clear that 


1 I ,2 P 

« pH + —7 

2 7 — 1 


= —div ( ^p\u\ 2 v^j -j-div(Prt) + u ■ divT. (4.6) 


From the continuity equation (11.11) i. momentum equations (II.ID ?, relation (|4.6D . Lemma 
m and integration by parts, we have 


—/(t) ——M(t ) — 2 (t + 1 )—F(t) + 2 (t + l) 2 —e(t) - 2 F(t) + 4 (t + l)e(i) 
(2 — 3(7 — l))(t + 1 ) f Pdx + Ji + J 2 , 


(4.7) 


7-1 

where Ji and J 2 are given by 


lA(t) 


Ji = —2(t + l) / x • divTdx, J 2 = 2(t + 1)^ 
JA(t) 


u ■ divTdx. 


'A(t) 


Since 


div(x • T) = x • (divT) + ^T a = x ■ (divT) + 3^/3 + p S divu. (4.8) 

Integrating (|4.8|) by parts, with help of the fact T = 0 on dA(t), we have 

J\ = — 2(t + l) f x ■ divTdx = 6(f + 1) f f (3 + f-a ) //divudx. (4-9) 
JA(t) JA(t) \ 3 J 
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Now we turn to Jo. From (11.61) and Cauchy’s inequality we have 


div(uT) =u ■ divT + 2 p(p) ^^(djitj ) 2 + p{p) ^^(<9, 




2=1 




+ 2 p{p) y ^(djUj)(djUi) + A (p) ( ^2 ®i u i 
i>j \i=l 

>u ■ divT + ^X(p) + (divu ) 2 

=u • divT + (/3 + J //(dm/) 2 . 


We integrate (14.101) over A(t) to get 


from which we have 


lA(t) 


j 2 < -2(t + iy 


' A(t) 


u ■ divTdx < — / ( /3 + -a ) //(dmt) 2 dx, 


•m 


(3 + -a ) p d (divu) 2 dx. 

O 


From Lemma ITTl along with (14.71) . (14.91) and (|4.12j) . for 0 < t < T m , we get 

d_ 

d t 

2 


f /(0<- ? - r (2-3(7-l))(t + l) / 
*7-1 

— 2(t + l) 2 


Pdx 


A(t) 

St 


f3 + -a ) p d (divu) 2 dx 

A(t) \ 3 


+ 6(f + 1) / ( (3 + -cm //dfvudx. 

4A(t) V 3 / 


(4.10) 


(4.11) 


(4.12) 


(4.13) 


A(t) 

Since 5 < 7 , with the help of Cauchy’s inequality and Young’s inequality, we have 
— 2 (t + l ) 2 / //(dfvu) 2 dx + 6 (t + 1 ) / //divudx 


’A(t) 


' A(t) 


< — 2 (i + l ) 2 / //(divu) 2 dx + 2(t + l ) 2 / p 5 (divu) 2 dx + 18 / //dx (4.14) 
JA(t ) 4A(t) J A(t) 


•A(t ) 

<18 [ p 5 dx < 18< ^ 


/A(t) 7 4A(t) 

We deduce from (14.131) that 

d_ 
dt 


[ p 1 dx+ 
JA(t) 


18(7 — 5) 


7 


A(t) 

|A)|- 


I(t) <—?—(2 - 3(7 - l))(t + 1) [ Pdx 
7 - 1 JA(t) 

+ 18 p 7 dx + 18 |A)|- 


(4.15) 
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From the second expression of I(t ) in (|4.5|> . one has 
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t + 1 


—— [ \x — (t + l)u\ 2 pdx 
t + 1 JA(t) 

+ -^-r(2-3(7-l))(t + l) [ 

7 - 1 JA{t) 


(4.16) 


Pdx. 


In the case when 1 < 7 < |, from fl4.15ll - fl4.16p . for 0 < t < T m , we have 


d_ 

dt 


2-3(7-1) 


<5(7 - 1) 


—I(i) <-7—j— -7(i) + 18 ( /3 + — 0 . ) -7-7— — —-— ^I(t) + 18 ( /3 + —a ) - |-A|- 


2A 1 {t+lf 


t + 1 

Solving (14.171) directly, we get 

J(t) <(t + l) 2 - 3(7 - 1) e-F 1 r ^ e Ol /(0) + a 2 J\t + i) 3 (7- 1 )- 2 e ^i d rj 


2 \ 7- (5 

A 1 

3 


7 


(4.17) 

(4.18) 


where 


«! = 18 | (9 + 


, 02 = 18 ( /? + -o^ ^ 


7 — (5 


Ail- 


ao(t + 1 ) _ai_ 

+ A e t +i e a i 


If 3(7 — 1) — 2 / —1, then from f[4.18p we get 

<C ^-A- 1 ) + t + 1) , for t € [0, T m ). 

If 3(7 — 1) — 2 = —1, from (14.181) we get 

I(t) <(t + i) 2 - 3 (7-i) e -^t ( e “i/(o) + a 2 e ai ln(t + 1)) 
((t + 1) ln(t + 1) + t + 1), for t <E [0, T m ). 


(4.19) 


(4.20) 


On the other hand, from the definition of I(t), Jensen’s inequality and Lemma 14.11 we 
show that 


I(t)> 2{t + ] )2 \A 0 \ f Ap\x,t) 
7 - 1 JA(t) 


dx 

'I^M 


> C(t + l )2 |A| 1 _ 7 m(0) 7 > C 0 (l + t) 2 , 
7-1 


(4.21) 


where Co > 0 is a constant and we used the fact in Lemma 14.31 that 

m(t) = / p(x,t) dx = / />o(x)dx = m(0). 

-/AO) -/An 


Then T m < +00 follows immediately, otherwise a contradiction forms between (14.211) and 
(ITT 7 )!) or Ml) . 
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In the case when | < 7 < + 00 , thus 2 — 3 ( 7 — 1) < 0* from (|4.15l) we have 

S /(t) S18 (/ + H t Lm P '' dX + 18 { l3+ I“) 


•m 

< 18 [p + ^a) + 


I(t) + 18 ( (5 + —o;^ l^oI 


ai 


r I(t) + a 2 . 


(4.22) 


(f + 1) 2 ' 

and therefore, 

I(t ) < e“ 1 /( 0 )e _ *+ I + a 2 e ai e _ ‘+Te ai t 
< e ai (7(0) + a 2 t). 

Again, this and (14.211) imply that T m < 00 . We complete the proof of Theorem 1.2. □ 


(4.23) 


4.2. Blow-up by hyperbolic singularity set. The mechanism for our second finite 
time blow-up result comes from the nonlinear hyperbolic structure (see (|1.151) 1 which 
controls the behavior of the velocity u in the vacuum region. Assume that the initial data 
( po,uo)(x ) have a hyperbolic singularity set V, see Definition 11.31 

Proof of Theorem 1.3. Let V(t) be the image of V under the flow map of (14.11) . i.e., 

V(t) = {x(t;{ 0 )\{ 0 eV}. (4.24) 

It follows from the continuity equation (ll.lD i that the density is simply transported 
along the particle path, so 

p(x,t) = 0, when x € V(t). 

From the Definition o for regular solutions, we have 

ut + u-Vu = 0, when x € V(t), (4.25) 

which means that u is a constant vector along the particle path x(t] £ 0 ) and 

£0 = x — tu{x, t ) € V. 

Then for any x € V(t), we have 

u(x, t ) = u 0 (£o) = u 0 (x- tu(x, t)), 


which implies that 

Vit(x, t) =(13 + t\7uo(x — tu(x, t))) 1 Vrio(x — tu(x, t)) 

= (I 3 + tVu 0 (Co)) _ 1 Vrto(Co), for x€V(t). 
According to the dehnition of the hyperbolic singularity set, there exists some 
£0 G V, and l ?0 € Sp(\7u 0 (£, 0 )) satisfying l$ 0 < 0. 

Let w € M 3 be the eigenvector of Vuo(£o) with respect to l^ 0 , that is, 

Vu 0 (£ 0 )w = l^w. 

It is clear that 

(l 3 + fVu 0 (£o)) -1 w = (1 + ti^ 0 y 1 w. 
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Thus we know that the matrix S7u(x,t) has an eigenvector w with the eigenvalue 

ho 

1 + tZf 0 ’ 

which, along with l^ 0 < 0, implies that the quantity Vtt will blow up in finite time, i.e., 

T m < Too. 

This completes the proof of Theorem 1.3. □ 


5. Appendix: Proof for the Remark 11.21 

In this section, we will show that the regular solution that we obtained in Theorem ll.il 
is indeed a classical one in (0,T*]. The following lemma will be used in our proof. 

Lemma 5.1. p.j If f(x,t) € L 2 ([0, T]; L 2 ), then there exists a sequence Sk such that 

Sk —>• 0, and Sk\f(x, s*;)|| — > 0, as k —> Too. 

From the definition of regular solution and the classical Sobolev embedding theorem, it 
is clear that 

(p, Vp, p t , u,Vu) € C(M 3 x [0, T*]). 

So it remains to prove that 

(ut, divT)(x, t) € C(M 3 x (0, T*]). 

From the proof of Theorem 1 1.1 1 in Section 3, we know that (through a change of variable 
s-i _ 

(p = P 2 ), system (11,11) can be written as 

4>t + u- V(p T ^ divu = 0, 

( \ (5-1) 

ut T u ■ Vu T -r- 4> 5-1 + (j> 2 Lu = V(j) 2 ■ Q(u). 

d — 1 

The solution (4>,u) satisfies the regularities in (13.61) and <f> € C 1 (M 3 x [0, T*])- 
Step 1: The continuity of ut■ We differentiate (15.11) 9 with respect to t to get 

u tt T 4> 2 Lu t = —((j) 2 )tLu - (u ■ Vu) t - T (V</> 2 • Q(u)) t , (5.2) 

which, along with (13.61) . implies that 


utt £ L 2 ([0,T*\;L 2 ). 


( 5 . 3 ) 
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Applying the operator di (|C| = 2) to (15.21) . multiplying the resulting equations by dim 
and integrating over R 3 , we have 


1 d_ 

2 dt 


\dim\ 2 + a\c/)V d^utll + (« + 


j ( ~ V< ^ 2 ' —J~ C ^ d x U t) - { d x(4> 

+ J did^tLu) -d£(u- Vu) 
+ J <9|(V^ 2 • Q(u)) t ■ dim 


+ / di(X7cj) 2 ■ Q(u)) t ■ dim 


- Vcj) 2 ■ 5 —^-Q{dim) - (di((j) 2 Lm) - 4> 2 Ldim)) ■ dim 


t — 



(5.4) 


15 



Now we analyze the terms J* (i = 10, • • • ,15). By Holder’s inequality, Lemma 12.11 and 
Young’s inequality, we have 



<q^V 3 ut|2|vV| 2 |V0|oo < C\m\ 2 D 2 + ^\<f>V 3 m\l 

■hi = J ~(di(c/) 2 Lut) - 4 > 2 Ldim) ■ dim 

<C^V 3 ^| 2 |V 2 ?xt|2|V0| oo + + |V 2 ^| 3 |(/)V 2 ^| 6 |r 4 | Z 52 ^ 

<C\m\ 2 D 2 + ^W 3 m\l 
■h 2 = J -di ((< i) 2 ) t Lu) ■ dim 

<C7 lool^tl^cz |<^V 4 Yi| 2 + |V 2 0| 3 |Lu| 6 |0 f | oo |R t | i: ,2 

+ |(/>V 2 'u t | 6 |(/)t| D 2|Lrt| 3 + |V^|oo|V(/)t|6|Lii| 3 |i{t| D 2 

+ |(/>V 3 rt| 6 |V(/>i| 3 |«t | D 2 + |V(/)| 00 |^| 00 |V 3 rt| 2 |'Ut| D 2 ^ 

<C(\Ut\ 2 D 2 + | 0 V 4 w |2 + 1 ^ + — | 0 V 3 Rt||, 

•h \3 = J ~di(u ■ Vu)t ■ dim 

<C'||u t || 2 ||u|| 3 |^| D 2 + / -(u-v)dim- dim 

^1 + \ut\j)2 + I V U \oo\dim\2^ < C^l + |'Ut|£)2^ , 
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J 14 = / - 


<c 


<9|v(y «-i • ajju* = J s -1 ^di -vdiut 


7-1 

2 i 


|V 2 ^| 2 |<£V 3 u t | 2 + |</>--%|</»t|oo|V^| 2 | 0 V ^| 2 
+ |0^^ 4 |oo|^|oo|V^| 6 |V^| 3 |^V 3 7x t | 2 + |^^- 3 |oo|V^ t | 2 |V^| oo |0V 3 ^ t | 2 ) 


(5.6) 


a 


<C + ^V\ t \l 


and 


Ji5 = J di(V(j) 2 ■ Q(u))t • d^u t 

<c(||V<)>|| 2 |u t |^2 + (||V0|| 2 |V«t| 3 + \\u\\ 3 U t \\ 2 )\0V 2 m\6 
+ (l|V^|| 2 |0V 3 « t | 2 + \\V^\\ 2 \cpV 2 ut\ 6 )\ut\ D 2 

+ (l|V^|| 2 ||0t|| 2 ||n|| 3 + ||^|| 2 |^V 3 n| 6 )|n t | Z) 2^ + J ^(V^ 2 ) t • Q(u ) • dim 
<c(\u t \ 2 D2 + |0V 4 n| 2 ^ + — |^>V 3 iq|| + ^151 ■> 

where 

■A51 = J ^x(V 0 2 )t • Q(u) ■ dim 

<C , ||V0|| 2 ||^|| 2 ||«|| 3 Kb2 + J WiVfr-Q(u)-diu t . 

Using integration by parts, the last term in (15.81) is estimated as: 

J 4>diV(i)t ■ Q{u ) • dim 

l^i|n 2 |V«|oo|«ib 2 + l0tb 2 |V 2 rt| 3 |(/>V 2 tt t | 6 + |^V 3 tt t | 2 |Vri| 


(5.7) 


(5.8) 


<C 


oo\m\\2 


(5.9) 


<C\u t \ 2 D 2 + —\(l)V 3 m\2- 

Then (15.41) reduces to 


o 3 lMl> 2 + —I^V 3 n t |2 < C(\u t \ 2 D 2 + |0V 4 n| 2 + 1) • 


, Ir , , lv .. „ IZ , (5.10) 

2 dt' ' u 2 lz ~ V ^ u J 

Multiplying both sides of (15.101) with s and integrating the resulting inequalities over [r, i] 
for any r E (0,7), we have 

+ J s\(j)V 3 m\2ds < CT\u t (T)\ 2 D 2 + C(1+ t). (5-11) 

According to the definition of the regular solution, we know that 

V 2 n t €L 2 ([0,T*];L 2 ). 

Using Lemma [5711 to V 2 ut, there exists a sequence Sk such that 

Sk —>• 0, and s^lV 2 tq(-, 'Sfc)| 2 —>• 0, as k —>• +oo. 
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Choosing r = Sk —» 0 in (15.111) . we have 

t\u t \ 2 D 2 + f s\cf)W 3 ut\lds < C(1+ t), 

Jo 

then 

t^Ut € L°°([0, T*]; iL 2 ). 

The classical Sobolev embedding theorem gives 

L°°([0, T]; H 1 ) n TT 1,2 ([0, T]; H^ 1 ) -a C([0,T];L 9 ), 
for any q € (3, 6]. From (15.31) and (15.131) we have 

tu t € C([0, T*]; VF 1 ’ 4 ), 

which implies that 

u t E C(M 3 x (0, T*]). 

Step 2: The continuity of divT. Denote N = 4> 2 Lu — S/<p 2 • Q(u). From equations (| 5.1 D o. 
regularities (13.61) and (15.13ft . it is easy to show that 

tN € L°°([0,T*];iL 2 ). 

Due to 

Nt E L 2 ([0,T*];L 2 ), 

we obtain from ()5. 141) that 

tN € C([0,T*]; IT 1 ’ 4 ), 

which implies that 

N E C(K 3 x (0,T*]). 

Since p E C(M 3 x [0, T*]) and divT = pH, we immediately obtain the desired conclusion. 

In summary, we have shown that the regular solution that we obtained is indeed a 
classical one in M 3 x [0, T*] to Cauchy problem (II.lft - (jl.3ft . 
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